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Abstract 

This paper studies countable systems of linearly and hierarchically interacting diffu- 
sions taking values in the positive quadrant. These systems arise in population dynamics 
for two types of individuals migrating between and interacting within colonies. Their large- 
scale space-time behavior can be studied by means of a renormalization program. This 
program, which has been carried out successfully in a number of other cases (mostly one- 
dimensional), is based on the construction and the analysis of a nonlinear renormalization 
transformation, acting on the diffusion function for the components of the system and con- 
necting the evolution of successive block averages on successive time scales. We identify a 
general class of diffusion functions on the positive quadrant for which this renormalization 
transformation is well-defined and, subject to a conjecture on its boundary behavior, can 
be iterated. Within certain subclasses, we identify the fixed points for the transformation 
and investigate their domains of attraction. These domains of attraction constitute the 
universality classes of the system under space-time scaling. 
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1 Introduction 



1.1 Model and background 

We are interested in the following system of coupled stochastic differential equations (SDE): 
dX v>i (t)= Yl ajv(£, rj) XvA*)] dt + ^J 2 9i{X v (t)) dB Vii (t), rj G Q N ,i = 1, 2. (1.1) 

Here a/v(-,-) is the transition rate kernel of a random walk on O/v, the hierarchical group 
(or lattice) of order N (see (jl.3p ). {X V } V ^ N with X v = (X V) i, X^^) is a family of diffusions 
taking values in [0, oo) 2 , g = (gi, g%) is a pair of diffusion functions on [0, oo) 2 , and {B TI } ve Q N 
with B v = {B„x, B n ^) is a family of independent standard Brownian motions on R 2 . As initial 
condition, we take 

X n (0) =0= (0 lt 2 ) G [0, oo) 2 Vry G U N . (1.2) 

Equation (jl.ip arises as the continuum limit of discrete models in population dynamics. In 
these models, individuals live in colonies labeled by the hierarchical group 0,^. Each colony 
r) G SI A? consists of two types of individuals, whose total masses are represented by the vector 
Xf,. Individuals migrate between colonies according to the migration kernel a/v(-, •). At each 
colony, each individual undergoes branching at a rate that depends on the total masses of the 
two types of individuals present at that colony. The system in (jl.lj) arises in the so-called 
"small-mass-fast-branching" limit, where the number of individuals in each colony tends to 
infinity, the mass of each individual tends to zero, and the effective branching rate grows 
proportionally to the number of individuals in each colony. The drift term in (jl.ip arises from 
the migration, which is the only source of interaction between colonies. The diffusion term in 
(jl.ip arises from the branching, where gi{x)/xi is the state-dependent branching rate of the 
i-th type, which incorporates the interaction between individuals within a colony. For more 
background, see e.g. Sawyer and Felsenstein |31j . Dawson and Perkins [16j . Chapters 9-10 in 
Ethier and Kurtz [21], Cox, Dawson and Greven [7], Dawson, Gorostiza and Wakolbinger [9]. 

The goal of the present paper is to study the universality classes of the large-scale space- 
time behavior of (jl.ip . It turns out that, for the specific form of the migration kernel a/v(-, •) 
given by (|1.5j) below and in the limit as JV -> oo, (jl.ip is susceptible to a renormalization 
analysis. The renormalization program for hierarchically interacting diffusions was introduced 
by Dawson and Greven [10], pj] for diffusions taking values in [0,1]. It has since been 
extended to several other state spaces (see Greven [26] for an overview). We will give more 
detailed references in Section \1. 31 First we outline the main ingredients of the renormalization 
program. 

1.2 Renormalization program 

The lattice in (jl.ip is the hierarchical group of order N, which is defined as 
n N = lr ] = (r ]l ) im £{0,l,...,N -1} N : ^r/i<cx)l 

I i£N J 

with coordinatewise addition modulo N. Define a shift <f> : O/y — > £In by (<ftr])i 
On fi/V) the hierarchical distance is defined as 

d(r], = mm{k G N = N U {0} : c/) k T] = (j) k ^}, 



(1.3) 

:= 7] i+1 (i G N). 

(1.4) 
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which is an ultrametric, i.e., d(r],£) < d(rjX) V c2(£, C) for all rj, £, C G ^jv- We choose the 
random walk transition rate kernel in such a way that a^(^, rf) depends only on the hierarchical 
distance between £ and r\. In view of what follows, we write in the form 



k>d(€,rj) 



N 



1-2* 



£,r? G ^at, £ 7^ 77, 



;i.5) 



where (c n ) ne N Q is a sequence of positive constants. Formula (|1.5p says that the random walk 
associated with ajv('> •) jumps with rate Cfc_i/iV fc_1 from 77 to an arbitrary site in the k-block 
{£ G fi n : </> fc £ = (j) k r]} around 77. 

The key objects in the renormalization analysis are the k-block averages: 



N k 



Q N , i = l,2, k e N . 



;i.6) 



Ceo, 



Using (jl.5p . we may rewrite (|l.ip as 



dX^Ct) = ]T Wp^{t) - X v>i (t)\dt + ^2 9i {X v {t)) dB Vji (t), r,en N ,i = 1,2, (L7) 



fc>i 



where each component X^ feels a drift towards the successive averages of /c-blocks containing 
n. It can be seen that the evolution of the 1-block averages is described in law by the SDE 



Ck 



fc>l 



— 9i(X^(tN)) dB Vii (t), rjen N ,i = l, 2, 



(1.8) 



where B v = {B r ,\,B„2} is a family of independent standard two-dimensional Brownian mo- 
tions. Note that in the limit N — > 00, we expect both the drift and the diffusion term in (|1.8p 
to be of order one, which means that evolves on the time scale tN. 

Let us next see heuristically what happens if we let N — » 00, the so-called hierarchical 
mean- field limit. If we let X — > 00 in (jl.7p . then the only drift term that survives is 

co[^£!i(*)-^,i(*)]^- 

Furthermore, Y^j(t) — > X(.)(0) = # for all i > 0, because evolves on the time scale tN. 

Therefore the system {X^i)}^^ converges in law to an independent system of diffusions, 
each satisfying the autonomous SDE 



dZi(t) = c (di - Zi) dt + J2 9i (Z(t)) dB.it) 



1,2. 



;i. 9 ) 



This kind of behavior is frequently referred to as "McKean-Vlasov limit" and "propagation of 
chaos" . 

With the above fact in mind, we move one step up in the hierarchy. Since Xg(t) evolves 
on the time scale t, for each fixed t the family 



{Xz(tN)}een N 



(1.10) 
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decouples and converges almost instantly to the equilibrium distribution of (jl.9p with the drift 
towards 9 replaced by a drift towards the first block average Y^\tN). Thus, we expect that 

Jj zZ 9i{X^{tN)) ~ J a T c ^ ]{tN) (dx) 9i {x) as iV -» oo for fixed t, 77 € f2jv, i = 1, 2, 

(1-11) 

where r?' 5 denotes the equilibrium distribution of ([1.90 . Thus, if we set 

(F co gW)= [ Tf 9 (dx) gi (x), i = 1,2, # € [0, oo) 2 , (1.12) 

then by (jl.lip . for large iV, the SDE (jl.8p for the 1-block averages Yjj takes exactly the same 
form as the SDE (jl.7p for the single components, provided that we rescale time by a factor 
N and replace the single component diffusion functions gi by (F Co g)i (i = 1,2). Here, F CQ 
plays the role of a renormalization transformation acting on the pair of diffusion functions 
9 = (51,52)- 

We can iterate the above procedure. The upshot of this is that, as N — > 00, the fc-block 
averages evolve on the time scale t N k according to the SDE 



dzf\t) = c k (Oi - Zf\t)) dt + ^2(^(11*]^)) dBi(t), i = 1,2, (1.13) 
with diffusion functions F^g = (F^gi,F^g2) given by 

^ = ^0-0^5, ken . (1.14) 

In fact, putting the successive iterates together and observing the sequence of block averages 

(YW(sN%Ylt? ] v (sN k ), ■ ■ ■ ,f}%N k )) (1.15) 

on the time scale sN k , as N — » 00, we expect this sequence to converge in distribution to a 
backward Markov chain 

(M(-k),M(-k + l),--- ,M(0)), (1.16) 
the so-called interaction chain, where 

(1) The starting position M{—k) is distributed as the weak solution of (|1.13p at time s with 
initial condition Z^(0) = 6; 

(2) for < j < k — 1, the transition probability kernel from M{—j — 1) to M(—j) is given 
by 

F[M(-j) G I M(-j - 1) = x] = Yf Fb]g [dy), (1.17) 

where 9 ( • ) denotes the equilibrium distribution of (|1.13p with k replaced by j. 

The distribution of M{— k) depends on s because Y^ (sN k ) evolves on the time scale sN k , 
while the transition probability kernel from M(—j — 1) to M{—j) for < j < k — 1 is 
independent of t because, conditioned on YYj +1 , YYj equilibrates almost instantly on the time 

scale sN k . Note that (F^g)i(9) = E[gi(M(0))\M(-k) = 9], where E denotes expectation 
with respect to the interaction chain. 

With these heuristics in mind, the renormalization program consists of the following two 
steps: 
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(I) Stochastic part: Show that for all scales fc E N, in the hierarchical mean-field limit 
N — > oo, the block average in (jf ,6p converges in law to the solution of the SDE in (|1.13p . 
and the sequence of block averages in (|1.15p converges in law to the interaction chain in 

(H3SD. 

(II) Analytic part: Analyze the renormalization transformation F c and the iterates F^ n \ 
n e N . 

Assuming that the stochastic part of the renormalization program can be completed, the 
large-scale space-time behavior of (II. ID in the limit N —* oo is characterized by the behavior 
of F^ as n — > oo, in particular, by its fixed shapes and their universality classes. 

Here, by a fixed shape we mean a pair of diffusion functions g = (51, 52) such that F c g = Xg 
for some c, A > 0. We speak of a downgoing fixed shape, fixed point or upgoing fixed shape 
depending on whether A<l,= l,or>l. Note that since the factor A can always be absorbed 
in time-scaling, such fixed shapes correspond to models that are mapped into themselves after 
a suitable rescaling of space and time. Indeed, if we set = c\ k (fc > 0), then such a fixed 
shape satisfies F^g = \ k g because the SDE associated with (cfc, F^g) is simply a time change 
of the SDE associated with (c,g), which induces the same renormalization transformation. 
For the interacting model in (|1.7p . this means that the fc-block averages evolve on the time 
scale tN k X k according to the diffusion function g. We note that our definition of a fixed shape 
deviates from the definition used in some earlier work, e.g. Fleischmann and Swart [23]. What 
is called a fixed shape there is, in our terminology, a joint fixed shape for all c > 0, i.e., a g 
such that for all c > there exists a A = A(c) with F c g = Xg. 

By a universality class, we mean a set Q of diffusion functions with the property that, given 
(cfc)fceN , for each g € Q there exist scaling constants (sn)neN such that s n F^g converges to 
the same limit (possibly up to a multiplicative constant). Typically, the limit will be a fixed 
shape or an asymptotic fixed shape (for the latter, see Fleischmann and Swart [24J). Note 
that each joint fixed shape gives rise to a universality class, namely all models within a given 
universality class exhibit the same large-scale space-time behavior. 

Apart from being relevant in the study of large-scale space-time behavior, fixed shapes 
also give rise to continuum models, by taking the so-called hierarchical mean-field continuum 
limit, which is a spatial continuum limit of the hierarchical lattice fijv with N — » 00. These 
continuum models also exhibit universality on small space-time scales, which is governed by 
the same renormalization transformation F c and its iterates F^ n \ n € No- For more details, 
see Cox, Dawson and Greven [7], and Dawson, Greven and Zahle |14j . 

The large-scale space-time behavior of (II. ID depends both on the diffusion function g 
and on the potential-theoretic properties of the random walk with transition rate kernel (jl.5p . 
Based on earlier work, we expect nontrivial universality classes to arise only when X^neN c n 1 = 
00, which is the "necessary and sufficient" condition for the random walk with transition rate 
kernel aj\r( - , ■) on fi^r to be recurrent (except for a side condition that becomes irrelevant 
in the limit N — > 00; see Sawyer and Felsenstein [31 J ) . For linear systems such as (jl.ip . 
the recurrence of the random walk is usually associated with clustering; see e.g. Dawson and 
Greven [IT], Cox and Greven [SJ, Swart [33]. In our context, clustering means that the solution 
of (II. ip converges in law to a mixture of distributions, each of which is concentrated on the 
configuration X-n = x, 77 € Qn, for some x 6 [0, oo) 2 with g\ (x) = g2{x) = 0. The choice 
of (cn)neNo determines the pattern of cluster formation, such as whether only small clusters 
appear, or only large clusters appear, or clusters of all scales appear. The latter is known as 
diffusive clustering (see e.g. Dawson and Greven |llj). 



5 



With the above facts in mind, the analytic part of the renormalization program can be 
more precisely formulated as follows. 

1 . Find classes of diffusion functions on which the renormalization transformations F c and 
their iterates F^ n \ n € No, are well-defined. 

2. Determine all the (asymptotic) fixed shapes. 

3. Determine the universality classes of diffusion functions that, for given {cn)n€N and 
after appropriate rescaling, converge to these (asymptotic) fixed shapes, and determine 
the associated scaling constants. 

1.3 Literature 

The full renormalization program has been successfully carried out for hierarchically interact- 
ing diffusions taking values in: 

(1) the compact interval [0, 1] (Dawson and Greven [10], [TT], Baillon, Clement, Greven and 
den Hollander [2]), where the Wright-Fisher diffusion is the unique fixed shape and is 
globally attracting with a scaling that is independent of the diffusion function; 

(2) the halfline [0, oo) (Dawson and Greven |12j . Baillon, Clement, Greven and den Hollan- 
der [3]), where the Feller branching diffusion is the unique fixed point and is globally 
attracting with a scaling that depends on the asymptotic behavior of the diffusion func- 
tion at infinity. 

For higher-dimensional diffusions, the analytic part has been carried out for: 

(3) isotropic diffusions taking values in a compact convex subset of R d (den Hollander and 
Swart [28], Swart [35]), where the diffusion function with constant curvature is the 
unique fixed shape and is globally attracting with a scaling that is independent of the 
diffusion function; 

(4) a class of probability-measure- valued diffusions (Dawson, Greven and Vaillancourt [13] . 
Dawson and March [15]), where the Fleming- Viot process is the unique fixed shape and 
is globally attracting with a scaling that is independent of the diffusion function; 

(5) a class of catalytic Wright-Fisher diffusions taking values in [0, l] 2 (Fleischmann and 
Swart [24] ) . where the diffusion function for the first component is an autonomous 
Wright-Fisher diffusion and the diffusion function for the second component is an au- 
tonomous Wright-Fisher diffusion function multiplied by a catalyzing function depending 
only on the first component. The renormalization transformation effectively acts on the 
catalyzing function. There are four attracting shapes for the catalyzing function, de- 
pending on whether the initial catalyzing function is zero or strictly positive at the 
boundary points of [0, 1] , and these attracting shapes are globally attracting with a 
scaling that is independent of the catalyzing function. 

The stochastic part for higher-dimensional diffusions has only been completed for interacting 
Fleming- Viot processes (Dawson, Greven and Vaillancourt [13]) and for mutually catalytic 
branching diffusions taking values in [0,oo) 2 (Cox, Dawson and Greven [7]). 
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All previous studies deal with diffusions that have certain simplifying properties. In the 
one-dimensional cases (1) and (2), as well as in the two-dimensional case (5), the equilibrium of 
(jl.9|) is reversible. As a result, many explicit calculations can be performed that are crucial for 
the analysis. For certain diffusions with compact state space, which includes the cases (1), (3) 
and (4) , there is a common underlying structure (called "invariant harmonics" , see Swart [34J ) 
that allows the determination of the unique fixed shape and its domain of attraction. In all 
cases where the state space is compact, the scaling needed for convergence to an attracting 
shape depends only on (c n ) ng N > n °t on the diffusion function g. This is different in case (2), 
where the state space is not compact. In all cases except case (5), the fixed shapes turn out 
to be joint fixed shapes for all c > 0. 

The goal of the present paper is to carry out the analytic part of the renormalization 
program for a general class of branching diffusions taking values in [0, oo) 2 . The multi- 
dimensionality and the non-compactness of the state space pose significant challenges. Due 
to the multidimensionality, the well-definedness of the renormalization transformation is non- 
trivial. The structure of the fixed points/shapes turns out to be rather rich. In fact, we 
will prove that, under certain restrictions, the class of fixed points is a ^-parameter family of 
diffusions with independent branching, catalytic branching and mutually catalytic branching as 
the extremal fixed points, and they are joint fixed points of F c for all c > 0. Moreover, we will 
prove that all diffusion functions that are comparable to these fixed points in an appropriate 
sense fall in their domains of attraction. 

1.4 Outline 

The rest of the paper is organized as follows. In Section [2] we formulate our main results, 
which come with varying degrees of restrictions on the diffusion functions. Section [3] contains 
the proof of the ergodicity of the SDE (jl.9p . and basic properties of the renormalization 
transformation. Section proves the identification of fixed points/shapes. Sections [5] and [6] 
identify the domains of attraction for the fixed points. In Appendices [A] and [B] we collect 
some technical results needed for the proofs. 

2 Main results 

In Section 12. 1\ we formulate a key class of diffusion functions C, for which the SDE (jl.9p 
has a unique weak solution. Section 12.21 contains a theorem on the ergodicity of the SDE 
(jl.9p . defines the renormalization transformation, formulates a subclass Tt + C C on which the 
renormalization transformation is well-defined and, subject to a conjecture on the preservation 
of certain boundary properties, can be iterated. Section 12.31 gives the definition of certain 
generalized fixed points /shapes, and identifies some special fixed points/shapes. Section l2~4l 
contains results on the identification of fixed points/shapes in Tt + under additional regularity 
assumptions. Section 12.51 contains our main result on the domains of attraction to the fixed 
points under further assumptions. Lastly, Section 12.61 provides a brief discussion of these 
results and lists some future challenges. 
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2.1 Key class and uniqueness for the autonomous SDE 



The renormalization transformation F c is based on (jl.9p . which is the SDE for the vector 
X(t) = (X 1 (t),X 2 (t)) G [0,oo) 2 written out as 



dX 1 (t) = c [0! - X^t)} dt + y /2g 1 (X 1 (t),X 2 (t))dB 1 (t), 
dX 2 (t) = c[6 2 - X 2 (t)\ dt + y/2g 2 (X 1 (t),X 2 (t)) dB 2 (t), 

where c > 0, 9 = {9\,9 2 ) € [0, oo) 2 , and B{t) = {B\(t) , B 2 {t)) are independent standard 
Brownian motions on M?. The corresponding generator is 

{LYfW)=cY j {e i -x i )-^J{x)+Y,9i{^m , /GC c 2 ([0,oo) 2 ). (2.2) 

i=l 1 i=l 1 

Note that, due to the absence of mixed partial derivatives, Ifl 9 can be interpreted as the 
generator of a two- type branching diffusion with state-dependent branching rates gi{x)/xi 
(i = l,2). 

Abbreviate 

Ai = [0, oo) x {0}, A 2 = {0} x [0, oo). (2.3) 
We will say that a function /: [0, oo) 2 — > [0, oo) has boundary property 



(d\) if lim = jiy) Vy € U A 2 with 7 continuous and > on A\ U A 



2. 



((?2) if li m _,4^ = j(y) Vy € U A 2 with 7 continuous and > on A\ U A2, (2.4) 



(^12) if lim^^^- = 7(2/) Vy G A\ U A2 with 7 continuous and > on A\ U A 2 . 

x— *y 1 2 

Throughout the paper, the pair g = (gi,g 2 ) will be assumed to be in the following class. 



Definition 2.1 [Class C] 

Let C be the class of functions g(x) = {g\{x) , g 2 {x)) satisfying: 

(i) For i = 1, 2, gi is continuous on [0, oo) 2 and > on (0, oo) 2 . 

(ii) For i = 1,2, gi satisfies boundary property (<9j) or {d\ 2 ). 



Note that for {g\,g 2 ) € C we can write gt{x) = Xi^/i(x) or gi(x) = x\x 2 ^i{x) for some positive 
continuous function 73 on [0, oo) 2 , depending on whether satisfies boundary property {df} or 
(S12). Note also that g\ and g 2 vanish on A 2 , respectively, A\, which is necessary to guarantee 
that the diffusion stays within [0, oo) 2 . Thus, if we denote the effective boundary of g by 

dg = {x£ [0,oo) 2 : 9l (x) = g 2 (x) = 0} , (2.5) 

then dg can be either of the following: 

Ax n A 2 , At, A 2 , AiU4 (2.6) 

These boundary constraints allow for the system (|2.ip to be treated as a perturbation of either 
of the following diffusions: 
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(1) Independent branching: (51,52) = (bixi, 02^2), 0i>&2 > 0, dg = Ai nyl2- 

(2) Catalytic branching: either (51,(72) = (oi^i, C2X1X2), 01, C2 > 0,dg = A2; or (51,52) = 
(01x1x2,62^2), ci,6 2 > 0,^5 = A\. 

(3) Mutually catalytic branching: (51,52) = (£1X1X2,02X1X2), ci,C2 > 0, dg = A\ U^2- 

Such a perturbation is behind the following result of Athreya, Barlow, Bass and Perkins p], 
and Bass and Perkins [B], which provides the starting point of our analysis. The latter paper 
developed out of Dawson and Perkins |17j . where Holder continuity is assumed rather than 
continuity, but the result there is not restricted to two dimensions as in [BJ. 

Theorem 2.2 [Well-posedness of martingale problem] (pQ, [6j) 

For all c> 0, 5 G C, 9 G [0, oo) 2 and x G [0, oo) 2 , with the possible exception of the case when 
x = (0,0), 9 G (0,oo) 2 , and either 51 or 52 satisfies boundary property {812), the martingale 
problem associated with the generator in (|2.2[) has a unique solution with starting position x. 

As a consequence of Theorem 12.21 the SDE (12.1D has a unique weak solution for all 9 G 
[0, oo) 2 and x G [0, oo) 2 , with the possible exception of the case when x = (0, 0), 9 G (0, oo) 2 , 
and either 51 or 52 satisfies boundary property (<9i2)- For each fixed 9 G [0, oo) 2 , the SDE 
(|2.1f) defines a Feller process satisfying the strong Markov property (see e.g. Theorem 4.4.2 in 
Ethier and Kurtz [21] and Corollary 11.1.5 in Stroock and Varadhan [33J ) . 

Remark 1: When 9 G (0,oo) 2 , 5 G C, 51 and 52 satisfy (<9i), resp. (82), the well-posedness 
of the martingale problem was established in Athreya, Barlow, Bass and Perkins pQ for all 
initial conditions x G [0, oo) 2 . When 9 G (0, oo) 2 , 5 G C, and 51, 52 both satisfy (<9i2), the well- 
posedness is established in Bass and Perkins [6j for all intial condition x G [0, oo) 2 \{(0, 0)}. 
Both [1] and [6J use local perturbation arguments and the results are not restricted to linear 
drift as considered here. Since the perturbation arguments are local, this implies that well- 
posedness also holds for mixed boundaries, i.e., 51 satisfies (d±) and 52 satisfies (^12), or vice 
versa. When either 51 or 52 satisfies (S12), Lemma 35 of Dawson and Perkins [17] shows that, 
for all x G [0, oo) 2 \(0, 0), with probability 1 the unique weak solution of ()2.ip with initial 
condition x never hits (0,0), and hence we can restrict the state space to [0, oo) 2 \{(0, 0)}. 
When 9 G 9[0,oo) 2 , the local analysis of pQ and [6j still applies until the diffusion first hits 
the absorbing boundary, at which time the diffusion becomes one-dimensional, a situation for 
which the well-posedness of the martingale problem is standard. 

Remark 2: The proof given in [lj requires the drift to be strictly positive in each component 
on 9[0,oo) 2 . However, as pointed out in Bass and Perkins [5], it is sufficient that the inward 
normal component of the drift is strictly positive on d[0, oo) 2 , which holds in our setting when 
0G(O,oo) 2 . 

Remark 3: It would be considerably more difficult to deduce from Theorem 12.21 the well- 
posedness of the martingale problem for the system (jl.ip . for which one would need to restrict 
the state space. To deduce the Feller property, one would need to restrict the state space even 
further and impose growth conditions on the diffusion function 5, typically 51 (x) + 52 (x) = 
0(x 2 +X2) (see, e.g., Shiga and Shimizu [32], Cox, Dawson and Greven [?])• We will not resolve 
these issues here, since they belong to the stochastic part of the renormalization program, 
which remains open. 
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2.2 Equilibrium distribution and renormalization transformation 

Our first result shows that (|2.ip has a unique equilibrium for the class C. The proof will be 
given in Section [3,11 Henceforth C denotes law. 

Theorem 2.3 [Equilibrium distribution] 

For all g G C, 9 G [0, oo) 2 and c > 0, (|2.ip has a unique equilibrium distribution T^ 9 ; which is 
continuous in 6 with respect to weak convergence of probability measures, and 

C(X(t)) =>■ T c f VX(0) G [0, oo) 2 . (2.7) 

The convergence in (12. 7j) is crucial for the stochastic part of the renormalization program (not 
considered here), while the uniqueness of the equilibrium is crucial for the definition of the 
renormalization transformation, which we now define. 



Definition 2.4 [Renormalization transformation] 

The renormalization transformation F c , acting on g G C, is defined as 

{F c g)i0) = [ 9i {x) Tf{dx), 6£ [0, oo) 2 , c> 0, i = 1, 2. (2.8) 

J[0,oo) 2 y 

Henceforth we will denote expectation with respect to T c l 9 by E^ 9 . 

6 

Without restrictions on the growth of g at infinity, it is possible that F c g is infinite. We 
therefore need to consider a tempered subclass of C. 



Definition 2.5 [Class Ho+] 

(i) For a > 0, let 7i a C C be the class of all g G C satisfying 

gi(x 1 ,x 2 ) + gi(xi, x 2 ) < C(l + xi)(l + x 2 ) +a{x\ + x%), (xi,x 2 ) G [0, oo) 2 , (2.9) 

for some < C = C(g) < oo. 

(ii) Let 

n 0+ = p| Ha. (2.10) 

a>0 



Note that Hq+ is much larger than Ho- In particular, Hq+ includes diffusion functions that 
along the axes grow faster than linear but slower than quadratic. 

Our second result shows that F c is well-defined on the class H a when < a < c, pre- 
serves the effective boundary, and preserves the growth bound in (|2.9p though with a different 
coefficient. The proof will be given in Section f3.21 



Theorem 2.6 [Finiteness, continuity, preservation of dg and growth bound] 

For c > and < a < c, if g G H a , then F c g is finite and continuous on [0, oo) 2 , dF c g = dg, 
and F c g satisfies (12. 9h with a replaced by ^S^a. 



To proceed with our analysis, we need the following: 
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Conjecture 2.7 [Preservation of boundary properties] 

Let g <E TC + . 

(i) For i = 1,2, if g^ satisfies then so does (F c g)i for all c > 0. 

(ii) For i = 1,2, if gi satisfies (812), then so does (F c g)i for all c > 0. 

In Section [3.31 we will explain why this conjecture is plausible. Combining Theorem 12,61 with 
Conjecture 12.71 we get: 

Corollary 2.8 [Preservation of class H +] 

For all c > 0, the class TCq+ is preserved under F c , i.e., F c g 6 TLq+ for all g G Hq+. 

The latter is a key property, because it allows us to iterate F c on T~Lq + and investigate the 
orbit F^ n 'g = F Cn _ 1 o - • ■ oF co g, n € No- We will not need Conjecture 12.71 or Corollary 12.81 until 
we study the iterates in Section 12.51 

The subquadratic growth bound imposed by T~Cq+ cannot be relaxed: we will see in Corol- 
lary 12.111 below that F c cannot be iterated indefinitely on TL a for any a > 0. 

2.3 Definition and examples of fixed points and fixed shapes 

We next give the definition of fixed points and fixed shapes of F c . Generalizing our definition 
given in the introduction, we allow for the case where F c g = \g with A not a constant but 
a diagonal matrix. These generalized fixed shapes do not give rise to universality classes as 
defined in Section 11.21 but they may be relevant for studying finer properties of the orbit 
(FW ff ) neNo . 

Definition 2.9 [Generalized fixed shapes and points] 

The pair g = (51,52) £ W a with a € [0, c) is called a generalized fixed shape of F c if 

Fc(9l,92) = (A151, A252) for some Ai,A 2 > 0. (2.11) 

If Ai = A2, then g is called a fixed shape, and if Ai = A2 = 1, then g is called a fixed point of 
F c . 

Our third result identifies a family of fixed points and (generalized) fixed shapes of F c . 
The proof is nontrivial because of integrability issues, and will be given in Section 13.21 

Theorem 2.10 [Examples of fixed points and fixed shapes] 

(i) The pair 

(91,92) = (hxi + C1X1X2, 62^2 + c 2 xix 2 ) (2.12) 

is a fixed point of F c in TCq+ for all c > and all b±, b 2 , ci, c 2 > with (61 + c\)(b 2 + c 2 ) > 0. 

(ii) The pair 

(91,92) = (a\x\ + bix\ + c\X\x 2 , a 2 x\ + b 2 x 2 + c 2 xix 2 ) ( 2 -13) 

is a generalized fixed shape of F c in 7^ ai va 2 f or all c > 0, < a\, a 2 < c and b\,b 2 ,c\,c 2 > 0. 
The corresponding scaling constants are 

Ai = — — , A 2 = — — . (2.14) 

c — a\ c — a 2 
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Diffusion functions of the form in (|2.12p are mixtures of independent branching, catalytic 
branching and mutually catalytic branching (recall Section 12. ip , all of which are in the class 
T~io+- We will see in Theorem 12,151 below that, under additional regularity conditions, such 
mixtures are the only fixed points of F c . Diffusion functions of the form in f)2 . 13|) are mixtures 
of these fixed points and the Anderson branching diffusion (51,52) = (oiscf, o^l)- The latter 
do not fall in the class Wo+- 

The following corollary of Theorem 12.101 shows that F c g cannot be defined for all g G TL a 
with a > c, and F c cannot be iterated indefinitely on 7i a for any a > 0. The proof will be 
given in Section 13.21 

Corollary 2.11 [Divergence of iterated fixed shapes] 

Let gi(x) = aixf+ fiiXi+^iX\X2 with 014 > and /3j,7i > 0, i = 1, 2. Let (c n ) n( z^ be the positive 
sequence that defines F^ n \ (see fli.i^D ). Let uq = min{n G N : (a\ Va^) ^27=0 c I l — 

((fW 5 ) 1 ,(fW« 7 ) 2 ) = (- ^r-r5i, ^r-T5 2 ), < n < n , (2.15) 

ictt (F[ n °]g)i + (F^g) 2 = 00 on (0,oo) 2 . 

2.4 Identification of fixed points and fixed shapes 

Our fourth result rules out generalized fixed shapes in Tt + with an upgoing component. The 
proof will be given in Section 14.31 



Theorem 2.12 [No fixed shapes in TLq+ with an upgoing component] 

For c > 0, there is no g G TCq+ such that either (F c g)\ = \\g\ with X\ > 1 or (F c g)2 = A252 
with A 2 > 1. 

Our fifth result does the same for generalized fixed shapes with a downgoing component, 
but only under mild additional regularity conditions. The proof will be given in Section [4.31 
Below, in line with general topological notation, liminf^_ > ( 00)00 ') denotes the infimum of all 
limits along sequences tending to (00, 00). 

Theorem 2.13 [Sufficient conditions for no downgoing fixed shapes in Hq+] 
Let c > 0. 

(i) There is no g G 7Yo+ such that F c (gi,g2) = (Ai<7i, A252) with < Ai,A2 < 1 and 

liminf (^ + ^)=0. (2.16) 

2^(00,00) \ xf x 2 J 

(ii) There is no g G Tlo+ such that (F c g)\ = \\g\ for some < Ai < 1 and g satisfies any of 
the following conditions: 

9l > on Ax \ {(0,0)}, (2.17) 
• liminf > 0. (2.18) 

x^(oo,oo) X\X2 

A similar result holds with the indices 1 and 2 interchanged. 
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Remark: Conditions (|2.16|) and (|2.18p are complementary. Note that one particular case 
not covered by conditions (|2.16f l (|2.18|) is when g\ vanishes on both axes, g± (x) = o(x\X2) as 
x —s- (oo, oo), and g 2 {x) = x\x 2 . In that case we cannot rule out the possibility of g\ being a 
downgoing fixed shape. 



In Theorem 12.101 we identified a 4-parameter family of fixed points. To show that these 
are the only fixed points, we need to impose strong additional regularity conditions. 

Abbreviate 

iioo = {(0,oo),(oo,0), (00,00)} (2.19) 

and 

fyoo.O) (^) = Xl > h (0,oa)(x) = x 2, fyoo.oo) (») = X \X 2 . (2.20) 

Definition 2.14 [Class H r ] 

Let TLq be the set of g G Tto satisfying 

(i) inf gi(x)>0 Vs>0, i = l,2, (2.21) 

xG[s,oo) 2 

(ii) lim = Xi g G [0, 00) Vz G Roo, i = 1,2. (2.22) 

hg(x) 

Note that Hq C Ho C 7io+- Also note that, because g\ vanishes on A 2 and g 2 on A%, 
necessarily A 1)(0)DO ) = A 2i(oOi0 ) = 0. 

Our sixth result is the following. The proof will be given in Section 14.11 

Theorem 2.15 [Identification of fixed points in Hq] 

LetoO and g = (gi,g 2 ) G U r Q . If F c (g 1 ,g 2 ) = {gi,g 2 ), then 



(2.23) 



gi (x) = A 1;(oo>0 )Xi + A 1)(o0)00 )aJiaj 2 , 

g 2 (x) = A 2i ( ,oo)2;2 + A 2 , (00,00) %l%2, 

where X^g, z£ i?oo, o,re defined in \2.22\ . 

2.5 Domain of attraction of fixed points 

Our seventh and final result is on the domain of attraction of the iterated maps 
F Cn _ 1 o ••• o F CQ , n G No, for a fixed positive sequence (c n ) ne N - We show that, provided 
inf ng N c n > and XmeN c « 1 = 00 > an diffusion functions that are comparable to a mixture 
of the fixed points fall into its domain of attraction. In Section [SJ we will give the proof for 
the special case c n = c, while in Section [6l we prove the result for varying c n . 

Theorem 2.16 [Domain of attraction of fixed points] 

Let (c n )neN Q be a sequence such that inf ng N c n > and X^neN c n X = 00 • Let g G H r be such 
that 

9i(x) > aiXi + PiXix 2 , ai,Pi > 0, + > 0, i = 1, 2. (2.24) 

Then 

lim (FWg)i(&) = V Xi,ghg($) V0 G [0,oo) 2 , i = 1,2, (2.25) 

n— >oo ' 

ZG-Roo 

where hg, Xig, z G -Rqo, are defined in \2. 20\) and \2.22\ . 
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What this says is that under the iterates F^ n \ any g that is properly minorized and has the 
same behavior at infinity as a mixture of the fixed points, converges to that mixture pointwise 
as n — > oo. 

Remark 1: Note that Theorem 12.161 implicitly assumes Conjecture 12.71 To be formally 
correct, in Theorem 12. 161 we should replace Hq by the largest subclass of Hq that is preserved 
by F c for all c > 0. 

Remark 2: The condition inf n6 N c n > means that we partially exclude the regime of large 
clusters (see e.g. Dawson and Greven [H]). We do not believe this assumption to be essential. 
As long as X^neN c « 1 = 00 > i- e -> the- associated random walk on Qn with transition rate kernel 
aj\r(-, •) is recurrent, we expect there to be universality and the convergence in (|2.25p to hold. 

2.6 Discussion and future challenges 

The results in Sections I2.2H2.5I constitute a partial completion of the analytic part of the 
renormalization program outlined in Section 11.21 We have formulated l~io+ as the class on 
which the renormalization transformation is properly defined and, apart from Conjecture 12.71 
can be iterated. We have proved absence of upgoing fixed shapes in this class, and absence of 
downgoing fixed shapes under mild regularity conditions, given by (j2.16j) - (|2.18p . Furthermore, 
we have identified our 4-parameter family of fixed points in (|2. X2j) as the only fixed points in 
a subclass Hq of the smaller class Ho, given by the strong regularity conditions (|2,2ip - (|2.22p . 
Finally, we have found the domain of attraction of these fixed points in Hq supplemented with 
the lower bound (|2.24p . i.e., diffusion functions that are comparable to a mixture of the fixed 
shapes. There are several open problems remaining, the chief among which are: 

(1) Verify Conjecture 12.71 i.e., establish that the renormalization transformation can be 
iterated on Hq+ . 

(2) Remove assumptions f)2. 16j) — (|2. 18j) in the proof of the absence of downgoing fixed shapes 
in H + . 

(3) Show that the fixed points in fl2 . 12f) are the only fixed points in Hq+. In particular, 
remove assumption A2.22 j> and the bound g\{x) + g2{x) < C(l + xi)(l + X2) in Hq C Ho. 

(4) Strengthen (2) and (3) by determining whether it is actually true that the fixed shapes 
in (|2.13p are the only fixed shapes in C. 

(5) Study the orbit of (F^g) n ^ when the behavior of g at infinity is different from that of 
the fixed points. In that case we still expect convergence, but only after F^g is scaled 
with n in some appropriate manner. For diffusions on the halfline [0, 00), this study was 
successfully completed in Baillon, Clement, Greven and den Hollander [3], which raises 
some hope that it can be carried through on the quadrant as well. 

The questions we treated in this paper and the open problems we just mentioned have close 
connections to probabilistic potential theory of diffusions and Markov chains taking values in 
the quadrant. Our proofs strongly lean on the observation that the fixed points we build 
are mixtures of extremal universal harmonic functions of the interaction chains described in 
Section 11.21 The problem of finding all fixed points then requires identifying the universal 
Martin boundary of these Markov chains. The reader interested in this point of view can 
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find the necessary concepts in Pinsky [30]. Harmonic functions have played an important 
role in earlier studies of the analytic part of the renormalization program. In particular, the 
convergence proofs in the cases (1), (3) and (4) mentioned in Section [TTBI all depend on a 
special property of these models, called "invariant harmonics" (see Swart |34j). Case (2) uses 
moment equations combined with comparison arguments, while case (5) uses a representation 
in terms of a superprocess. Due to multi-dimensionality and non-compactness, these tools 
either do not apply or are insufficient for our model. However, our present methods have 
their limitations as well. In particular, in their present state they can only be used to prove 
convergence to joint fixed points of F c for all c > 0, as opposed to fixed shapes, or cases where 
there might be different fixed points of F c for different values of c. Moreover, we can treat 
only functions that are perturbations of these fixed points, albeit in a rather large class. 

Another interesting question is to study multi-type branching models with more than two 
types. The class of random catalytic networks introduced in Dawson and Perkins [17J and 
generalized in Kliem [29] provide a rich class of fixed points of the renormalization transfor- 
mation. However our results here do not extend trivially to higher dimensions, because we 
need the well-posedness of the martingale problem (Theorem I2.2|) . which is more delicate in 
higher dimensions. Also, our proof of the formula (|A.3|) for the mixed moment X1X2 does 
not extend to mixed moments of higher order. 



3 Proof of Theorems [2J2, \2Jg &7W\ and Corollary I2JH 

In Section 13.11 we give the proof of Theorem 12.31 in Section 13.21 that of Theorems 12.61 12.101 
and Corollary 12.111 In Section 13.31 we discuss Conjecture 12.71 Along the way we need a 
proposition on moment equations for the equilibrium distribution r^.' 9 , which will turn out to 
be fundamental in our analysis. This proposition is formulated and proved in Appendix lAl 



3.1 Proof of Theorem 12^1 

We break down the proof of Theorem l2.3l into four parts: existence, uniqueness, weak continuity 
and convergence. For uniqueness and convergence, we need to distinguish between 6 € (0, oo) 2 
and 6 € <9[0,oo) 2 . 
Existence: 

Proof. If we denote the distribution of X{t) by fit, with fiQ = 8$ for some arbitrary x G [0, oo) 2 , 
then it suffices to show that \yt : vt = \ Jq fi s ds}t>o forms a tight family of distributions on 
[0,oo) 2 . Indeed, we can then find a sequence (t n ) tending to infinity such that u tn converges 
weakly to a limiting distribution v. Consequently, for any / € C 2 ([0, oo) 2 ), 

(L c ff)(x)u(dx) = lim !{Lff)[S)u tn {dS) 

lim-?-/ I "{Uz 9 f)(x)n s (dx) ds 
(L c ff)(X(s))ds 



n^oo t n j 

lim -3-E M0 

n—*oo t, n 



= lim -5-E w [/(X(i n )) - f(X(0)} = 0, (3.1) 

n— >oo t n 

where the first line uses that Ut n converges weakly to v, the second line uses the definition of 
ut n , the third lines uses the definition of jJL s and Fubini, and the fourth line uses that f{X{t)) — 
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f(X(0)) — Jq(L c J 9 f)(X(s))ds is a martingale and / is bounded. Since J {Lfz 9 f)(x)u{dx) = 
for all / G C 2 ([0,oo) 2 ), which form an algebra of functions that is dense in the space of 
continuous functions on [0, oo) 2 vanishing at oo, it follows from Theorem 4.9.17 in Ethier and 
Kurtz [21] that v is an equilibrium distribution for (|2.ip . 

Tightness of the family {vt]t>o follows from the following lemma. 



Lemma 3.1 [Tightness estimate] 

Let (X(t))t>o be the unique solution of the martingale problem for L c ^ 9 with initial condition 
X(0) = x. Then 

E [X t {t) - 6i] < (xi - 6i)e- ct , i = 1, 2, t > 0. (3.2) 
Proof. For any pi, p 2 > 0, the function f(t, x) = Y^.=i Pi( x i ~ &i) eCt satisfies 

(Lf + -)/(*, x) = cY,Pi(0i - x t )e ct + cYjnlxi - 6^ = 0, (3.3) 

i=l i=l 

and therefore the process Yli=i Pi(Xi(t) — 6i)e ct is a local martingale. Introduce stopping 
times 

T n = inf \ t > : ^pjXj(t) > n I , n £ N. (3.4) 



Then 

2 



^Piixi-Oi) = Y;P^[( X ^ tAT n)-8)z c{tATn) ] (3-5) 

i=l i=l 

2 2 

= Y.P^K^-^hr^t}} +Y J P^[(Mr n ) ~ 9 l )e CT "l {Tn < t} ] . 
1=1 1=1 

For n > Ylf=i P$ii * ne second term in the right-hand side is nonnegative, so letting n — > oo 
we find that 

2 2 

J^PiE [X;(t) - f ] e ct < J2Pi( x * - 9 i)- ( 3 - 6 ) 

i=l i=l 

Since pi,P2 > are arbitrary, we arrive at (|3,2p . I 
This completes the proof of the existence. I 
Uniqueness: 

Proof. We distinguish between 6 in the interior resp. on the boundary of [0, oo) 2 . 

G (0,oo) 2 : By Theorem^ the unique weak solution (X(t))t>o of (|2.ip is a strong Markov 
process. By Remark 1 following Theorem l2.21 we restrict the state space to be [0,oo) 2 \{(0,0)} 
for the cases where weak uniqueness is not known when X(0) = (0,0). If (X(t))t>o has two 
distinct equilibrium distributions, then we can find two extremal equilibrium distributions p 
and v that are singular with respect to each other (see e.g. Theorem 6.9 in Varadhan |36j). 
This implies that there exist x, y G [0, oo) 2 such that the transition kernels Pt(x, •) and pt(y, ■) 
are mutually singular for all t > 0. However, if x, y G (0, oo) 2 , then we can first apply Theorem 
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IB, 41 to transport the diffusions started at x, resp. y, to a common small neighborhood with 
positive probability, and subsequently apply Corollary IB.3I to see that pt(x,-) and Pt(y,-) 
cannot be singular for all t > 0. On the other hand, when either x or y E <9[0, oo) 2 , it suffices 
to note that the drift in (|2.ip forces the diffusion to enter (0, oo) 2 instantly, which we justify 
shortly. Then, again by Theorem IB.41 the diffusion can be kept in (0, co) 2 up to any fixed 
time with positive probability, which reduces it to the case x, y E (0, oo) 2 . 

We now show that, for X(0) = x E <9[0, oo) 2 , (X(t)) t > enters (0,oo) 2 instantly. Consider 
first the case X(0) E {0} x (0, oo). Let X(0) = (0,y) with y > 0, and let r e = M{t > 
0: \X 2 (t)-X 2 (0)\ > y/2 or X x (t) > e}. Then Xi(tAr e ) -J^ T ' c(6q -X x {s))ds is a martingale, 



\2{L)—sl2{V)\ y/6 or ss*i\L) d. tf. men SL\\b/\T e ) — $ ' 

and 

" ptAT e 

(3.7) 



E[Xi(tAr £ )] =E 
Letting t — > oo, we find that for e small, 

e > E[Xi(r e )] =E 



/•tAr 6 

/ c(6i -Xi(s))ds 

JO 



3(6>i - Xi(s))ds 



> ^E[r e ]. (3.8) 



Therefore E[r e ] — > as e | 0, which is possible only if (X(t))t>o enters (0, oo) 2 instantly. The 
case X(0) £ (0, oo) x {0} is analogous. For X(0) = (0,0), a similar argument shows that 
(X(t))t>o enters [0, oo) 2 \{(0, 0)} instantly, which reduces it to the previous cases. 

9 E d[0, oo) 2 : If 6\ = 0, then E^ 9 ^!] = 0\ = for any equilibrium distribution Y c -1 9 by 

$ 6 

Proposition IA.ll In particular, T^' 5 is concentrated on {0} x [0, oo). Furthermore, (X\(t)) t >o 
is a local super martingale, and hence {0} x [0, oo) is an absorbing set. The equilibria for 
(X(t))t>o are therefore exactly the equilibria for (X(t))t>o restricted to the axis {0} x [0, oo), 
which is a one-dimensional diffusion. The proof of the existence and the uniqueness of the 
equilibrium distribution for this one-dimensional diffusion can be deduced either from explicit 
calculations as in Baillon, Clement, Greven and den Hollander [3J, or from the same argument 
as above for the two-dimensional diffusion with 9 E (0, oo) 2 . The situation is similar if 9 2 = 0. 

■ 

Weak continuity: 

Proof. We will show that T^ 9 is weakly continuous in 9. Let (9 n ) be a sequence such that 
9 n — > 9 in [0, oo) 2 . It suffices to show that {r c ; 9 } ng N is tight, and that any weak limit point of 
rl' 9 is an equilibrium distribution for the SDE (|2.1|) . which must be the unique T ! 9 . Tightness 
of {r^ ,5l } ne N follows from (|A.2|) . Suppose that T c ^ 9 converges weakly to a distribution v. Then 
for any / E C 2 ([0, oo) 2 ), 

(Lf f)(x)v{dx) = f {Ll 9 f){x)T c f(dx) + [ \(Lf - L c f)f\ (x)T c f{dx) 

[0,oo) 2 ^ J[0,oo) 2 y " °" J[0,oo) 2 L y y ™ J tin 



+ / (L c ff){x)[v{dx)-T c f{dx% (3.9) 



where the first term is zero because T c ^ 9 is an equilibrium distribution for the SDE in (|2.ip 

with parameter 9 n , the second term tends to as n — > oo because / E C 2 ([0, oo) 2 ) and 
\\L C ^ 9 f(x) — L C Z 9 f(x) ||oo — > as 9 n — > 9, and the third term tends to as n — > oo by the weak 

convergence of T c f to v. Therefore f, Q ^{L^ 9 f)(x)v(dx) = for all / E C 2 ([0, oo) 2 ). By 
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Theorem 4.9.17 in Ethier and Kurtz [21], it follows that v must be an equilibrium distribution 
for (|2.ip . and hence v = T c -1 9 . I 



Convergence: 

Proof. We again distinguish between 9 in the interior resp. on the boundary of [0, oo) 2 . 

9 E (0,oo) 2 : Firstly, note that by Theorem IB.4I and the fact that (X(t)) t >o started from 
<9[0, oo) 2 enters (0, oo) 2 instantly (see the paragraph containing (|3.7H3.8p above), the equilib- 
rium distribution T'h 9 must assign positive measure to every open subset of (0,oo) 2 . 

6 

Secondly, we show that for almost all x € [0, oo) 2 with respect to r^ 9 , C(X(t)\X(0) = x) 

9 

converges weakly to rl' 9 as t — > oo. We achieve this by showing that, for almost all (x, y) € 

6 

[0, oo) 2 x [0, oo) 2 with respect to the product measure T c 4 9 x T^' 9 , we can couple two solutions 

6 6 

(X(t)) t >o and (Y(t)) t >o of §2J]) starting from x, resp. y, such that lim*-^ F(X(t) ^ Y(t)) = 0. 
This goes as follows. 

Let e,S > be chosen as in Corollary IB. 31 where b(x) = c(9 — x) and a(x) = ^ 91 q^ ff2 (#)) 

on [0, oo) 2 (the definition of (a, b) in the rest of the plane M? is irrelevant, for instance one 
may define it by reflection), D = {x £ [0,oo) 2 : \\x — (1,1)|| < |} and x* = (1,1). Note 
that a(-) is nondegenerate on D for g € C. If (X(t))t>o, (Y(t))t>o are two independent 
copies of the strong Markov process defined by (|2.ip . then the joint process (X(t),Y(t))t>o 
is strong Markov and, by the same argument as for a single diffusion (X(t))t>o, the joint 
process has a unique equilibrium given by the product measure T- 9 x T- 9 , which implies that 

6 

the stationary process (X(t), Y(t)) t >o with £(X(0),Y(0)) = T C J 9 x T^ 9 is ergodic (see e.g. 

— o 

Theorem 6.9 in Varadhan [36] and the remarks thereafter). Since T- 9 x T^ 9 assigns positive 

6 6 

measure to B e (x*) x B e (x*), by the ergodic theorem almost surely (X(t),Y(t))t>o visits the 
set B e {x*) x B e (x*) after any finite time T. In particular, for almost all (x, y) with respect 
to T^ 9 x T^ 9 , almost surely the Markov process (X(t),Y(t))t>o starting from (x,y) visits 

u u — 

B e (x*) x B € (x*) after any finite time T. For such a pair (x, y), we construct the coupled process 
as follows. Start the independent processes (X(t)) t >o and (Y(t)) t >o with initial conditions x, 
resp. y. Thenr = inf{t > : (X(t),Y(t)) S B € (x*)xB e (x*)} < oo almost surely. By Corollary 
IB.3| the conditional transition probability kernels [ij* = F(X(t + 5) € • |(X(r), Y(r))) and 
= P(Y(t + 5) S ■ \(X(t), y(r))) have a common part H-^-p with measure at least \. 
From [i£ x jx^, we can take out fijff x /J,j£-&, which has measure at least | , and couple 
(X(t + 5 + t))t>o and (?(r + 5 + i))t>o so that they coincide for all t > and evolve as 
the strong Markov process defined by (|2.ip with initial measure l^vy- With respect to the 

remaining measure fXjt x ^ — /x^ ^ x ^ ^, we let (X(r + 5 + i))t>o and (Y(t + 5 + t))t>o 
continue to evolve independently. Since fijt x /i^ — Hfy x ^xy ^ s absolutely continuous with 

respect to fj,jt X /i^, a.s. (J? (r+<J+t)j (Y(r+(5+t))t>Q will visit B e (x*) x B e (x*) again. We can 
therefore iterate the above coupling procedure. Each iteration reduces the probability that X 
and Y have not been successfully coupled by a factor -j. Continue the iteration indefinitely to 
get the desired coupling between X and Y. We comment that, unlike in the context of Harris 
chains (see e.g. Section 5.6 of Durrett [18]) where one would need F(X(5) £ -\X(0) = x) to 
be dominated from below by a positive measure uniformly for x € B e (x*), to get a successful 
coupling it suffices that F(X(5) G -\X(0) = x) and F(X(S) € -\X(0) = y) overlap with 
probability at least a for some a > uniformly for all x, y £ B e (x*). 
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Next we show that, for Lebesgue almost every x G [0, oo) 2 , C(X(t)\X(0) = x) =^> T^ 9 

as t -> oo. Let A = {x G [0,oo) 2 : £(X(i)|X(0) = x) i=> T c f\. By Theorem O and the 
remark following it, the process defined by (12, ip is Feller continuous, and therefore ^4 is Borel- 
measurable. If A has positive Lebesgue measure, then we can find a simply connected bounded 
open domain D C (0, oo) 2 with smooth boundary such that AnD has positive Lebesgue mea- 
sure. We have shown above that T C -^ 9 (A) = 0, and hence T c -l 9 (A n D) = 0. If (X(t))t>o is 

the stationary solution of (|2.ip with marginal distribution T^ 9 , then E[j l^me^n-D^l = 
for all T > 0. On the other hand, by Theorem IB.51 we have for every x G -D that 
E[J rD lx(t)eAn£)^ I -^(^) = > 0. Since T c ^ g assigns positive probability to D, we have 



rt' 9 (dx) > o. 



By the monotone convergence theorem, we can choose T sufficiently large such that 

rr D AT 



E 

D 



X(0) 



r c /(dx) > o, 



the left-hand side of which is in turn dominated by E[f^ ^x(t)eAnD^ = ^> wmcn is a contra- 
diction. Therefore A has Lebesgue measure 0. 

Lastly, we show that C(X(t)\X(0) = x) T^ 9 for all x G [0,oo) 2 . Indeed, for x G 

9 

(0, oo) 2 , let e > be such that B e (x) C (0, oo) 2 . By Corollary HOI applied to D = B e (x), the 
transition kernel /if e ^\x,-) with killing at the boundary of B e (x) is absolutely continuous 
with respect to Lesbesgue measure. Since, for Lebesgue almost every y G B e (x), C(X(t + 
s)\X(t) = y) T^ 9 as s — > oo and [if^ x \x, B e (x)) | 1 as t | (see (IB.3P ). we have 
£(X(t))\X(Q) = x) => T^ 9 . The case x G d[0, oo) 2 follows from our previous observation 
that X(t) starting from x enters (0, oo) 2 instantly (see (|3.7p - (|3.8p ). 

9 G d[0, oo) 2 : Without loss of generality we may assume that 9\ = 0. If -Xi(O) = 0, then 
Xi(t) = for all t > and (Xi(i), A^(i)) = (0, A^i)) is effectively a one-dimensional diffusion 
with diffusion function 32(0, x%). By the same argument as before, albeit much simpler, this 
one-dimensional diffusion is ergodic, and the convergence in (|2,7p holds. If -Xi(O) 7^ 0, then it 
suffices to show that Xi(t) — > a.s. and £(A^(i)) =>- T^ 9 as i — ► 00, where T^ 9 is taken as a 
measure on [0, 00). 

Note that Xi(t) is a local supermartingale and X\(t) A 1 is a bounded supermartingale, so 
that X\{t) A 1 — > y a.s. as i — ► 00 for some non-negative random variable Y. By the bounded 
convergence theorem and (|3.2p . 

E[Y] = lim EpTi(i) A 1] < lim Xi(0)e~ c * = 0. (3.10) 

t^QO t—>00 

Therefore Y = and Xi(t) —* a.s. as f — > 00. 

To show that C(X2(t)) ==> T c 4 9 as t ^ 00, it suffices to show that E[<^>(A^2(i))] — > 

6 

E c -i 9 [4>(X2)] as t — > 00 for any </> G C 2 [0, 00). Abbreviate 

6 

a = E c J 9 [<f>(X 2 )] and u(t, x) = E{<P(X 2 (t)) \ X(0) = x\. (3.11) 

9 

For X(0) G [0,oo) 2 with Xi(0) = 0, (X(t)) t > is effectively a one-dimensional diffusion that 
is ergodic, and hence u(t, x) — ► a as t — > 00 for each x G {0} x [0, 00). We claim that in fact 
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u(t, x) — > a uniformly on compact intervals of the form {0} x [0, if]. To see why, note that if 
Y{t) and Z{t) are solutions of the one-dimensional SDE 



dX{t) = c{9 2 - X(t)) dt + y / 2g 2 (0,X(t)) dB t (3.12) 

with initial condition Y(0) = y < Z(0) = z, then Z{t) stochastically dominates Y(t) for 
all t > 0, i.e., if F ttV (v) = P(Y(t) < v\Y(0) = y), then F %y (v) > F ttZ (v) for all t,v > 0. Let 
Foo(v) = T c J g (— oo, v). Then, for any x 2 > 0, F t X2 (v) — > -Foo(^) as t — > oo for all but countably 
many v G [0, oo). For any x 2 G [0, if], if > 0, we can write 



«(t,(0 J x 2 ))=/ <p{v)dF tjX2 (v) = - <f>'(v)F t ^(v)dv= (cf>'_(v)-(f>+(v))F t>X2 (v)dv, (3.13) 
io io JO 

where <j/ + (v) = <j/(v) V and <f>L(v) = -(</)' (y) A 0). Since 

i f'OO f'OO 

<t>L(v)F t)K (v)dv < / <j>l(v)F t7X2 (v)dv < / <fL(v)F t , (v)dv, (3.14) 



where both ends of the inequality tend to (/)'_(v)Foo(v)dv by the bounded convergence 
theorem, <j)'_(v)F tjX2 (v)dv converges uniformly to J °° ^^Fooiv^dv for x 2 G [0, if ] as t — > 
oo. A similar statement holds for J Q (f>'+(v)Ft tX2 (v)dv. Therefore u(t,x) converges uniformly 
to a on {0} x [0,K\. 

Let X(Q) G [0, oo) 2 be arbitrary. By (|3.2p . (AT 2 (i))j>o is tight, and hence for any e > we 
can choose if large enough so that ¥(X 2 (t) > if) < e for all i > 0. Since u(t, x) — > a uniformly 
on {0} x [0, if], we can choose ti large enough so that supa. 2e [o,j<] l u (*ii (0, x 2 )) — a| < e/2. 
Since { (X(i)) t >o}^ ( )e[0 oo) 2 defines a Feller process (see the remark below Theorem 12. 2p . 
u(t\,x) is continuous in x G [0,oo) 2 . We can therefore choose 5 > sufficiently small so that 
sup^ g [ 0i5 ] x j 0i _R-] \u(tx,x) — a\ < e. Since Xi(i) — > a.s., we can choose t 2 large enough so that 
¥(Xi(t) > 5) < e for all t > t 2 . Then, by the Markov property, for any t > t\ + 1 2 we have 

u(t,X(0))=E[u(t 1 ,X(t-t 1 ))} 

w(ti,^(<-*i))l^ (t _ tl)6[ o ia]x[0 ,K] (3.15) 



E 



+ E 



u(ti,X(t i l))lx( t _ tl )^[o i «5] x [0,X] 



Since P(X(t — ii) ^ [0,(5] x [0,if]) < 2e and |M|oo < Halloo, a < IHloo, we easily verify from 
(I3TT51) that 

\u(t, X (0)) -a\ < e + 4e||0|[oo for all t > h + i 2 - (3.16) 
Since e > is arbitrary, u(t,X(0)) — > a as t — > 00, and hence £(X(t)) =^> T^ 9 . I 



3.2 Proof of Theorems EH f2TT0l and Corollary I2TTT1 

Proof of Theorem EH Let g = (gi,g 2 ) G H a for some < a < c. Then, by (pll . there 
exists a < C = C(i?) <oo such that 

9i{x) +92{x) < C7(l + ari)(l + x 2 ) +o(zf + a^), (xi, x 2 ) G [0, oo) 2 . (3.17) 

The finiteness of -F c <7 follows from Proposition lA.lf ii). If # ra — > for some G [0, oo) 2 , then, 
by Proposition lA.ll iii). <7i,g 2 are uniformly integrable with respect to {r^ s } rag N. Combining 
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this with the fact, shown in Theorem 1 2 . 3 1 and proved in Section \3. 11 that T c -l 9 converges weakly 

On 

to Yf as 9 n -> 6, we have E c /[ gi {X)) -> E^[ 5i (X)], i.e., (F c g)i(9 n ) -» (F c g)i(9) for i = 1,2 
(recall (|2~%|) ). 

By the moment equations (|A.2HA.3j) . we have 

(F^)i(e) + (i^) 2 (0) = E^bi(X) + 52 (X)] 

< E c f[C(l + X x )(l + X 2 ) + a(X 2 + X 2 )] (3.i 8 ) 

= C(l + 6>i)(l + a ) + a(# 2 + el) + ^((F c5 )i(0) + {F c g) 2 {9)). 

Therefore 

(F^e) + (f c9 ) 2 ($) < -^(c(i + 9 1 )(i + e 2 ) + a(e 2 1 + e 2 2 )). (3.19) 

c — a V / 

Consequently, if g G W +, then F c <7 satisfies (|3. 19[) for all a > 0, and so it satisfies the 
subquadratic growth bound imposed by the class H + . 

To show dF c g = dg, note that F c g > is obvious. If 9 € (0, oo) 2 , then the equilibrium 
distribution T c ^ 9 has positive mass in (0, oo) 2 , and so (F c g)(9) > follows from the fact that 
g > on (0, oo) 2 . If 9\ = 0, then, by (|A.2[) . T C J 9 is concentrated on the vertical axis A 2 . Since 

gi vanishes on A 2 , it follows that (F c g)\(9) = 0. Moreover, {F c g) 2 {9) = if and only if g 2 
vanishes on A 2 (recall (j2.5H2.6fl ) . A similar result holds for 9 2 = 0. I 

Proof of Theorem 12.101 Theorem 12 . 101 (i) follows immediately from (|A.2HA.3]l . To prove 
Theorem ETU](ii), note that, by (|X2lfA~4|) . 

{F c g) x {9) = E c J 9 [ ai Xf + hXt + Cl XxX 2 ] = ai E c f[X 2 } + Mi + c x 9 x 9 2 

= a 1 9 2 1 + -(F c g) 1 + b 1 9 1 + c 1 9 1 9 2 = g 1 (9) + ^(F c g) 1 (9). (3.20) 
c c 

Solving for (F c g) 1 (9), we get (F c g)i(9) = -^gi(9). Similarly, we have {F c g) 2 = -^g 2 for 
g 2 = a 2 x\ + b 2 x 2 + c 2 xix 2 . The assumption (bi + c\){b 2 + c 2 ) > is meant to rule out the 
uninteresting case g\ = or g 2 = 0. I 



Proof of Corollary 12. Ill Equation (|2.15j) follows from Theorem l2.10l (ii) by induction. Note 
that if a« X^fclo 1 c fe 1 — ^ or e ither i = 1 or 2, then the coefficient of x\ in (F^ n °~^g)i(x) is 
«</[! - "iEfcio 2 ^ 1 ] ^ c no-i- To show (F 1 ™ ^)! + {F [no] gh = oo on (0,oo) 2 , it therefore 
suffices to show {F c g) \ + (F c g) 2 = oo on (0, oo) 2 for g of the form gi{x) = a{X 2 + fiiXi +^iX\x 2 
with a.\ V a 2 > c. Without loss of generality, assume a± > c. The proof of Proposition lA.ll (ii) 
shows that the moment equations (|A.2HA.4j) are valid as long as {F c g)i{9) + (F c g) 2 (9) = 
E c ^ 9 [gi + g 2 ] < oo. Assume (F c g)i(9) + {F c g) 2 (9) < oo for some 9 G (0, oo) 2 . Then 

E^[X 2 ] = 9\ + -(F c5 )!(0) = 0? + ^E c f[X\\ + ^ + ^9 X 9 2 , (3.21) 
e c c c c 

which is not possible for oi\ > c. Therefore we must have (F c g)i(9) + (F c g) 2 (9) = oo for all 
0e(O,oo) 2 . I 
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3.3 Discussion of Conjecture 12.71 



In this section we explain why Conjecture 12.71 is plausible. We focus on the case where gi,g 2 
both satisfy boundary property (di 2 ) in (12. 4ft . i.e., gi{x) = xix 2 r yi(x) and g 2 {x) = x\x 2 ^ 2 {x) 
with 71,72 > continuous on [0, oo) 2 . 

Consider the tilted equilibrium 

t^m = j^ T 7( d ^ ^ (°^) 2 ' ( 3 - 22 ) 

where ()A.3[) implies the proper normalization. The conjecture amounts to showing that, as 
9 — > 9* G <9[0, oo) 2 , this tilted equilibrium converges weakly to some probability distribution 
on [0, oo) 2 , say T c J 9 (dx), that is weakly continuous in 9* and, in addition, 7 i(x) is uniformly 

integrable with respect to T*2 9 {dx) for 9 in a small neighborhood of 9*. Indeed, this observation 

6 

is immediate from the identity 

f ll (x)t^(dx) = -±-(F c g) l (9), i = l,2. (3.23) 

Now, recalling the generator in (|2.2p . we note that V'z 9 {dx) is the equilibrium associated 

9 

with the time-changed diffusion given by the generator 



c{e x -xi) d c(9 2 -x 2 ) d d 2 d 2 

(Lff)(x) = — / (x) + — f(x) + 7l (i)— -/(irJ + TaCa;)-— 2/(0?) 

/ £ C 2 ([0, oo) 2 ), 0-x)- V/(x) = on 8[0, oo) 2 . (3.24) 



Let 9^9* = (a, 0) for some a > 0. Then, at least heuristically, we get a limiting 
generator 

- - + + (3 . 25) 

/eCV 2 ([0,oo) 2 ). (9*-x)-V/(f)=0 on 8[0, oo) 2 \{9*}, JLf(S>) = JLf(S') = a. 

Here, the diffusion part has no singularity at the boundary, but the drift part does. As the 
process approaches the vertical axis A 2 it feels a growing drift downwards and to the right, 
while as it approaches the horizontal axis A\ it feels a growing drift horizontally towards 
(a, 0) and a constant drift downwards. Therefore, again heuristically, this generator describes 
a process that is obliquely reflected in the direction of (a, 0) upon hitting A 2 , and upon hitting 
A\ jumps to (a, 0) instantly and then moves back into the interior by reflection. Like the 
original diffusion with generator (I2.2p . this process ought to exist, be weakly unique, and have 
an ergodic equilibrium T C J? that is weakly continuous in 9* G [0,oo) 2 . 



4 Proof of Theorems [27T2L [27131 and [27151 

Section T4. II contains the proof of Theorem 12. 151 which is an immediate consequence of Propo- 
sition 14.11 Section 14.21 contains some preliminary lemmas needed for the proof of Proposition 
14.11 Section 14.31 provides the proof of Proposition 14.11 and of Theorems 12.121 and 12.131 
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4. 1 Proof of Theorem I2TT51 



The proof of Theorem 12.151 is based on an asymptotic analysis of the homogeneous Markov 
chain M c,g = (M c ' g (n)) n ^ with transition probability kernel given by p(9,dy) = T c -! 9 (dy), 

the unique equilibrium distribution of (|2.1|) . For F c g = g, M c,g is in fact the interaction chain 
in (|1.16p . Throughout the rest of the section, unless specified otherwise, we will denote the 
Markov chain M c,g by X. For F c g = g and g G both gi and gi are harmonic functions of 
X, i.e., both (gi(X (n))) n ^ and (g2(X (n))) n ^ are martingales. Theorem 12. 151 then follows 
immediately from the following proposition. 

Proposition 4.1 [Harmonic functions of X = M c,g ] 

If g G H.Q+ and satisfies ()2.21[) in the definition of Hq, then every nonnegative harmonic 
function f of X = M c ' 9 , i.e., every f such that 

E[/(X(n)) | X(0) = 6] = f{9) V0 G [0,oo) 2 , n G N , (4.1) 

which furthermore satisfies the constraints 

(i) f{x) < C(l + si)(l + x 2 ) /or some < C = C(/) < oo, (4.2) 

(ii) lim/(x) = Vzedg, (4.3) 



(m) = A /i? G [0,oo) VzGiioo, (4.4) 



Ijm ; urn 
zs of the form 

f(x) = ^2 X f,z h z(%) = A /,(oo,0) x l + "V,(0,oo) ^2 + ^,(00,00) XiX 2 , (4.5) 

ZS-Roo 

with hg, zeRoo, given by ([2TT9]) (f2T20|) . 

The proof of Proposition l4.1l will be given in Section l4~3l The strategy is to first /i-transform 
X (see Definition 14.31 below) to a new process X h = (X h (n)) n& fq using 

h(x) = (l + Xl )(l+x 2 ), (4.6) 

i.e., X h is defined as the homogeneous Markov chain with transition probability kernel 

p(e,dy) = h(y)T c J 9 (dy)/h(9), 

which is well-defined since h{x) is a harmonic function of M c,g . The function / is harmonic 
for M c ' 9 if and only if f /h is harmonic for X h . The constraint in (|4.2j) guarantees that f jh 
is bounded, the constraints in (|4.3H4,4p guarantee that f /his continuous up to the boundary 

R = 8gUR 00 , (4.7) 

while the constraint in (|2.2ip guarantees that lmin.^oo X h (n) G R a.s. It is then standard to 
show that f jh is uniquely determined by its values at R, which will imply (14.50 . 

The proofs of Theorems 12.121 and 12.131 are also based on an asymptotic analysis of the 
Markov chain M c ' g , even though when g is not a fixed point of F c , it no longer corresponds 
to the interaction chain in (|1.16p . 
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4.2 Preliminary lemmas 

The key results in this section are Proposition 14.61 and Corollary 14.71 

Let X = M°' 9 be as stated before Proposition 14. 11 First we list some moment equations 
for X(n), n £ No, which follow immediately from Proposition I A . 1 [ 

Lemma 4.2 [Moment equations for X = M c > 9 ] 

Let c > 0, and g G 7i a for some < a < c. Fix X(0) = 9 € [0, oo) 2 . Then for all n G N , 

E[X i {n)]=9 i , i = l,2, (4.8) 

E[X 1 {n)X 2 (n)} = 9 1 9 2 . (4.9) 

// ((F c g)i, {F c g) 2 ) = {Xigi, X292) for some Ai, A 2 > 0, then 

ng i (X(n))]=\ 1 ?g i (9), i = 1,2, (4.10) 

1 n 

nxf(n)} =9* + -J2 X hi(9), i = 1, 2. (4.11) 

In the proof of Theorem 12.15^ we will need Doob's /i-transform of a Markov chain, which 
we recall here. For more information on the /i-transform, see e.g. Section 4.1 of Pinsky |30j . 

Definition 4.3 [/i-transform] 

Let X = (X(n)) n£ N be a Markov chain with state space E and n-step transition probability 
kernel p n (x,dy). If h is a nonnegative (not identically zero) harmonic function of X, i.e., 
(h(X(n))) n £N is a nonnegative martingale, then the h-transform of X, denoted by X h , is 
defined as the Markov chain on the space {x G E: h(x) > 0} with n-step transition probability 
kernel pfo(x, dy) = p n (x, dy)h(y) / h(x) . 

The next two lemmas are immediate consequences of Definition 14.31 
Lemma 4.4 [Harmonic functions of X h ] 

Let X, h and X h be as in Definition \A.3[ If f is a harmonic function of X, then f/h restricted 
to {x G E: h(x) > 0} is a harmonic function of X h . The converse is true if h(x) > for all 

xeE. 

Lemma 4.5 [Absolute continuity of X h w.r.t. X at bounded stopping times] 

Let X, h and X h be as in Definition IP! If X(0) = X h (0) = x G E where h(x) > 0, and r 
is a bounded stopping time, then the law of X h (r) is absolutely continuous with respect to the 
law of X(t) with density ^4. 

The next proposition is the key to establishing Proposition 14. 11 Such a result is referred to as 
almost sure extinction versus unbounded growth, see e.g. Fleischmann and Swart |24j . 

Proposition 4.6 [Almost sure limit of /i-transform of X = M c ' 9 } 

Let c > 0, and let g G Hq+ satisfy condition ()2.2ip . Let h(x) = (1 + xi)(l +X2) and let X h be 
the h-transform of X . Then, for any X h (0) G [0, oo) 2 , almost surely, lim n ^oo X h (n) = X h (oo) 
exists and X h (oo) G R (see fl^.Tfl ). 
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Before giving the proof of Proposition 14.61 which we defer to the end of this subsection, we 
first state and prove a corollary and another prerequisite lemma. 



Corollary 4.7 [Trapping probabilities] 

Let c, g, h, X h and X h (oo) be as in Proposition ^.^ 

OO 

^(0)^(0) 

:i +A?(o))(i + *£(())) 

(ii) If (0, oo) x {0} ^ dg, then 



P[**(oo) = (oo, oo)] = - ; ^ZnlUr^ ' ^ 



P[^ h (oo) = (oo,0)] = j— . 



(4.13) 



(Hi) If {0} x (0, oo) ^ dg, then 



P[X h {oo) = (0, oo)] = rMr^ r--T. (4.14) 

L V ' V n (1 + ^(0))(1+X£(0)) V ; 

Proof By Lemmas 14.21 and 14. 4[ 

t l ~*\ XlX2 t (~\ Xl f (~\ X2 (A 1 C\ 

(1+Xl)(l+X 2 ) (1+Xl)(l + X 2 ) (1 + Xl)(l+X 2 ) 

are bounded harmonic functions of .X" , and therefore (/j(X h (n))) ne N , i = 1,2, 3, are bounded 
martingales. Since, by Proposition 14.61 X (n) — > X^(oo) E R a.s. asn-> oo, we have 

fi(X h (0)) = E c *[f i (X h (w))], i = 1,2,3. (4.16) 

Now (|4.121f47l"4"j) follow from the following observations: (1) /i((oo,oo)) = 1 and f\ = on 
i?\{(oo,oo)}; (2) if (0,oo) x {0} £ dg, then / 2 ((oo,0)) = 1 and f 2 = on i?\{(oo,0)}; (3) if 
{0} x (0, oo) £ dg, then / 3 ((0, oo)) = 1 and f 3 = on R\{(0, oo)}. I 

The proof of Proposition 14.61 in turn relies on the next lemma, which gives a lower bound 
for t : 9 (dx) = V2 9 (dx)h(x)/h(9), the transition kernel of X h with h(x) = (1 + xx)(l + £ 2 ), 

that is uniform in both g and 9. The uniformity in g is not needed for the proof of Proposition 
6J but will be crucial for the proof of Theorem 12.161 in Section [5j 



Lemma 4.8 [Uniform lower bound on T C J 9 (dx)] 

Let A C H + . 

(i) For any 9 € [0, oo) 2 , if 



3 e' > such that inf gAx) > for i = 1 or i = 2 (4.17) 

9S-4 

ses e ,(9) 

2 . M^f 



with B e ,(0) = {xG [0,oo) 2 : ||:r-0|| < e'}, then 



3 e > such that inf Vh 9 h ( [0, oo) 2 \B e (9) ) > 0. (4.18) 

2es e (e) 



25 



(ii) For any a > 0, if 



3 e', N' > such that 



inf #2 (x) > 

x£[N' ,oo) X [(* — e^a+e'] 



(4.19) 



and 



Va > 0, 3 C a € [0, oo) suc/i £/iai, uniformly for all x € [0, oo) 2 and g £ A, 
gi(x) + g 2 (x) < C a (l + xi)(l + x 2 ) + a(x\ + x 2 ,), 



(4.20) 



then 



Be, N > such that inf r5 s h ([0, oo) 2 \[iV, oo) x [a - e, a + e]) > 0. (4.21) 

x£ [JV,oo) X [a — e,a+e] 



A statement similar to f \4-21 ) holds for vertical strips of the form [a — e, a + e] x [N, oo) if, in 
l \4- 19\ , g 2 is replaced by g\ and [N r , oo) x [a — e' , a + e'] is replaced by [a — e', a + e'] x [iV 7 , oo). 



Proof. We first prove (|4,18p and (|4.2ip with replaced by T^' 9 . The main tool is the 
following moment equation valid for g £ 7Yq+> # £ [0, oo) 2 and i = 1, 2: 



E' 



1 



(1 + ^) 2 



1 + a 



■E ! 



<-:() 



1 



1 + Xi 



+ 



c(l + 9i) 



E" 



(1 + X<) 



(4.22) 



where X = (X(t))t>o i n this proof denotes the stationary solution of the SDE (12. ID . By 
stationarity, £(X(s)) = T c -1 9 for all s > 0. Hence 



Mi(t) 



rt L cJ 1 

i+Xi(t) i+Xi(o) y ?v!+^ 



ds, i = 1,2, 



(4.23) 



are local martingales, where 



9 







x=X(s) 



^2 



9 



L^ 9 = 0(0! - + c ( 2 -x 2 )^- + 0l (f)^-2 + <7 2 (x)^|. (4.24) 

Since E c /[Xi(s)] = 9 { and E c f [ gi (X (s))] = {F c g)i(6) < oo by Proposition EH we have 



E 



sup \Mi(s)\ 

L 0<s<* 



< 2 + E 



(--,9 



L J0 



c\9 i -X i (s)\ + 2g i (X(s)))ds 



< 2 + 2t(c6 i + (F c g) i {6)) < oo. 



(4.25) 



Therefore Mi = (Mi(t))t>o, i = 1)2, are in fact martingales, and E^ s [Mj(i)] = 0. By the 
stationarity of X, we have 



E 



(-41 



L 



<--<j 



1 



1 + Xi 



E 



eg 



1 + Oi-l-Xi 2 9i (X) 
-c • - — h 



x=X(s) 

Rearranging terms, we obtain (|4.22p . 



(1 + Xi) 2 (1 + X % 



0, i=l,2. (4.26) 
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(|4.18p : Suppose that (|4. 18j) with T^ 9 h replaced by T^ 9 is false. Then 

inf r c J 9 (\0,oo) 2 \BJe)) = Ve>0. (4.27) 

By (|4,17p . we may assume without loss of generality that mi g& ^^ £B ,^,g\{x) = 5 > for 
some eo > 0. In particular, mi gG ^g &B ,^gi(x) > 5 for all e G [0, eo]. Fix e G [0,eo]. Let 

G B e 0) and £ A be chosen such that ([0, oo) 2 \B t (9)) = o(l) as n -> oo. In 

P~22j) with i = 1, substitute x^ and #( n ) for 6 and 5. Then 

Lh - s --(iT^F + o(1) ' , , 

1 2 5 ^ 

r - h - s ' - (ITWT^) + WTeT) x (1 " 0(1)) x (i + 6> 1 + e )3 ' 

where we applied Jensen's inequality to obtain Tjq^TJ 2 * n ^ ne es ti m &te for the r.h.s. For 
e > sufficiently small and n sufficiently large, the above two equations are incompatible, and 
therefore (|4TT8|) with t c 4 g h replaced by T c f holds. Since h(x) = (l + x 1 )(l + x 2 ) > 1 on [0, oo) 2 

and is bounded on B t (9), it is easy to see by the definition of t^H that (|4. 18[) also holds . 

P~2Tj) : The proof that (|43Tj) holds with T c / h replaced by T C J 9 is the same as above and we 
leave the details to the reader. To get (|4.21j) . we argue as follows. 

Choose e G (0, a) and Nq > such that 

a , e , No = inf T c f ([0, oo) 2 \[7Vo, 00) x [a - e, a + e]) > 0. (4.29) 

x£[Nq,oc) X [a — e,a + e] 

By Proposition IA.1[ we have 

E^[(Xi -X!) 2 ] = ^-(F.gUx), g G H 0+ , x G [0, oo) 2 . (4.30) 

Therefore 

F c /{y G [0, oo) 2 : m < xi/2} < F c J 9 {y G [0, oo) 2 : \ Vl - xi\ > xt/2} < ^h^) (431) 

CX -y 

We claim that 

lim sup V , V ' = 0. (4.32) 

Assume ()4.32p for the moment. Since f3 a ,e,N is nondecreasing in N, we can choose > A^o 
sufficiently large such that 

inf r^{y*G[0,oo) 2 \[iV,oo)x[a- e ,a + e ]: y x > ^} > ^f^. (4.33) 

[iV,oo) X [a — e,ai + e] 
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Then 



inf f£J ([0, oc) 2 \[iV, oo) x [a - e, a + e]) 

r£ [iV,oo) X [a — e,CK + e] 



inf 



mNM-e, a+ e] ^[0,oo)2\[JV,oo)x[a~ e , a +e] 



> 



inf 



)(i +y-2)„ , 



9 eA J( !/l>*lA \ (1 + Xl)(l + X 2 ) 

S6[iV,oo)x[Q-e,a + e] \ [0,oo) 2 \[JV, oo) X [a-e,a+e] / 



> 



inf 



l + xi/2 



> 



(l+Xl)(l+Q + e) 7/ Wl>*l/2. 1 

S6[iV,oo)x[a- f ,c« + f ] 1 [0,oo)2\[JV,oo)x[a-e,Q + e] / 

4(l + a + e) ' 



which establishes ()4.2ip . 

To verify (|4.32p . note that, by condition (14.201) and Proposition IA.H 

E^bi + g 2 \ < [C a (l + Xt)(l + X 2 ) + a(X? + X|; 

= C B (1 + xi)(l + x 2 ) + a{x\ + x 2 ) + -E<f + g 2 ] V g € A. 



(4.34) 



(4.35) 



Solving for Eg[gi + g 2 ], we get 



E c J 9 [ gi +g 2 } = (F c g) x {x) + {F c g) 2 {x) < (c a (l + x x )(l + x 2 ) + a(xj + 



Therefore 



4(F c #)i(x) 4ca 
limsup sup g — • 

X\ — >oo g€A CX-^ C (X 



Since a > can be made arbitrarily small, (|4.32p follows. 



V g E A. 
(4.36) 



Proof of Proposition 14. 61 By Lemma 14.21 h\{x) = 1 + x\, h 2 {x) = 1 + x 2 and h(x) = 
(1 + xi)(l + x 2 ) are harmonic for X. Hence, by Lemma 14.41 h\(x)/h(x) = 1/(1 + x 2 ) 
and h 2 (x)/h(x) = 1/(1 + xi) are harmonic for X h . Therefore (1/(1 + X^(n))) ng N and 
(1/(1 + X^ - (n))) nG pj are nonnegative martingales and, by the martingale convergence theo- 
rem, X (n) — > X ft (oo) G [0, oo] 2 a.s. as n — > oo. We need to show that 



(i) P X h {oo) G [0,oo) 2 ,X /l (oo) £ d 5 



0. 



Xf(oo) = oo,X£(oo) G (0,oo) =P X^(oo) =oo,X^(oo) G (0, oo) 



0. 



If (i) is false, then there exists a # G [0,oo) 2 \dg such that, for all B C [0,oo) 2 with 
9 G int(S), F[X h (n) G B for all n large enough] > 0. In particular, we must have 



inf 



x&B h(x) 7[o i00 )2\B 



h(y)T c f(dy) =0 VBc[0, oo) 2 with 9 G int(B). 



(4.37) 



Otherwise, there is a uniform probability of escaping from B at each step, and X h cannot be 
confined in B forever with positive probability. 
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If (ii) is false, then (considering without loss of generality the first part of (ii)) there exists 
an a G (0, oo) such that 



P X h (n) G [N, oo) x [a - e, a + e] for all n large enough > Ve G (0, a), N > 0. (4.38) 
In particular, we must have 

^ inf -L / h{y)Tf{d$) = Ve > 0, JV > 0. (4.39) 

xe[N,oo)x[a-t,a+t] tl{X) J[o,oo) 2 \[AT,oo) x [a-e,a+e] 

But both (I4.37P and (|4.39p contradict Lemma [4.8l applied to A = {g}, where conditions (|4.19l - 
I4.20p in Lemma 14.81 are easily verified by our assumption that g G 7i + and that g satisfies 
(12.21j) . Therefore we must have lim n _ >00 X h (n) = X h (oo) G R a.s. I 



4.3 Proof of Proposition 14.11 and Theorems 12.121 and 12.131 

Proof of Proposition |4?0 Let / be a nonnegative harmonic function of X = M c,g satisfying 
the constraints in (|4.2H4,4p . Since x\, x 2 and x\X 2 are harmonic for X, so is fo(x) = f(x) — 
^f,(o,oo) x 2 — ^/,(oo,o) x i ~~ ^f,(oo,oo) x i x 2- Let X h denote the ^-transform of X with h(x) = 
(1 + x\)(l + x 2 ). Then, by Lemma 14.41 fo/h is harmonic for X h , and so 



fo{9) _ E t fo(X h (n)) 
h{6) L h{X h (n)) 



X h (0) = 9 V n G N, 6 G [0,oo) 2 . 



(4.40) 



Constraint (j4.2[) implies that fo/h is bounded, constraint (|4.4p implies that lim^^ fo(x) j h(x) 
= for all z G Rqq, while constraints (|4.3H4.4p imply that lim^g fo(x)/h(x) = for all 
z G dg. S incG, by Proposition 14. 6[ lini^^oc, X h (n) = X h (oo) G R(= dg U R oo) a.s., letting 
n — > oo in (|4.40p and applying the bounded convergence theorem, we obtain fo/h = and 
/o = 0. Therefore f(x) = Xf^o,oo) x 2 + ^f,(oo,o) x l + A/, (oo,oo) x l x 2- ' 

Proof of Theorem I2.12t Suppose the claim is false. Then, without loss of generality, we 
may assume that ((F c g)i, (F c g) 2 ) = (Ai5i,A252) for some g G W +, Ai > 1, Ai > A2 > 0. 
By Definition 12.51 f° r an y a > there exists a < C a < 00 such that g\(x) + g 2 (x) < 
C a (l + xi)(l + x 2 ) + a(x\ + x\). Fix X(0) = 6 E [0, oo) 2 , then by LemmaO we have 



E[ 5l (X(n))] 



< 



E 



C«(l + Xi(n))(l + X 2 (n)) + a(X 2 (n) + X 2 2 (n)) 



(4.41) 



< C a (l + 0i)(l + 2 ) + a(fl 2 + 2 2 ) + - E ( A i5i(^*) + A2W)) • 

C j=i 

Since Ai > 1 and Ai > A2 > 0, dividing both sides of the above inequality by A" and letting 

n — ► 00, we get 



9i(9) < 



c(Ai - 1) 



9i{d) + lx 1= x 2 92(9) 



(4.42) 



Since a > can be made arbitrarily small, (|4.42|) implies that g\(6) < 0, which is a contra- 
diction. I 
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Proof of Theorem 12. 13t (i) Assume that, for some g G Tt + with liminf^( 00 00 )[g , i(x)/a;j + 
g 2 {x))/x 2 2 ) = 0, F c ( gi ,g 2 ) = (Ai 5 i, A 252 ) for some < Ai, A 2 < 1. Fix A(0) = G [0,cx)) 2 . By 
Lemma l4.2[ we have 



1 n \ 
E[(AQ(n) - Oi) 2 } = -J2 A^(0) < - . <k(0), 



Vn G N. 



(4.43) 



Next, choose such that e ^' 1 _ A 4 < ^ for i = 1, 2, which is possible by the above assumptions. 
Then, by the Chebychev inequality, 

X(n) G 



01 30i 

~2~'~2~ 



02 302 

~2' T~ 



and hence 



E 



1 

> - 



2 ^r e i 38 



inf 



e 2 3e 2 



X(0) = 

> 



> 



Vn G N, 



Vn G N, 



(4.44) 
(4.45) 



which contradicts the assumption that K[gi(X(n))] = \fgi(0) — > as n — > oo. 
(ii) We consider the conditions (|2.17p and (|2.18|) separately. 

(|2.17p : Assume that (F c g)i = Xigi with Ai < 1 and gi(x\,0) > for all x\ > for some 
g G TCo+- For # = (#i,0) with #i > 0, T c -l 9 (dx) is supported on the horizontal axis A\ and is 



in fact the equilibrium distribution of the one-dimensional diffusion 



dXi(t) = c(0! - Xi)dt + y /2g 1 {Xi,0)dB 1 (t). 



(4.46) 



Therefore the mapping g\(xi,0) i— > (F c g)i {x\, 0) is the renormalization transformation for 
diffusions on the halfline which, by Lemma 2 and Theorem 2 in Baillon, Clement, Greven and 
den Hollander |3j, cannot have a fixed shape with scaling constant Ai ^ 1. 

t|2. 18j) : Assume that {F c g)\ = Xigi with Ai G (0, 1) for some g G W + such that lim inf^->(c« oo) 
g\{x)/x\X2 = e > 0. Then the /i-transformed Markov chain X h with /i(x) = (1 + x\){\ + x 2 ) 
satisfies E[( 5 i//i)(A fe (n))] = X r l(g l /h){X h (0)). If A^(0), X£(0) > 0, then, by Corollary Ml 

Xf(0)X£(0) 



PpT(oo) = (oo,oo)] 



(1 + A^(0))(l+At(0)) 



> 



and 



lim A" 



> 



h(xHo)) 

X?(0)X 2 h (0) 



lim E 

n— >oo 



9i(X h (n)) 
h{X h {n)) 



lim inf > , 



(1 + Zf (0))(1 + (0)) 5-*(oo,oc) /i(x 



(4.47) 

(4.48) 
(4.49) 



which is a contradiction. 
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5 Proof of Theorem 12.161 with constant 



Proof. Assume c n = c > 0, in which case 

F [n] = F n The proof is based 

on an analysis of the 

interaction chain introduced in Section H. 21 Let g satisfy the conditions in Theorem 12 .161 Let 
X = (X(—n)) n £^ be the (inhomogeneous) backward Markov chain on [0, oo) 2 with transition 
probability kernel 

X(-n) G dx X(-n - 1) = 6] = T c -l F " 9 (dx). (5.1) 



Denote the transition probability kernel from time — m to time — n > —m by K~ m ~ n (x, du). 
By Proposition IA.ll the functions 1, x\, xi and x\x 2 are harmonic for X. Let X h = 
(X h (—n)) ne ^ denote the /i-transform of X with h(x) = (1 + x\)(l + x 2 )- Then 1, ^ , 



x 2 



1+X 2 



— and T7^v are harmonic for X h . Now change variables and let 

X2 h[x) ° 



Y(-n) = cP(X h (-n)), 



with 4>: [0, oo) 2 — > [0, l) 2 given by 



</>{xi,X2) 



Xi 



x 2 



1 + x 1 1 + x 2 



(5.2) 



(5.3) 



Then Y = (Y(— n)) n eN i s a backward Markov chain on [0, l) 2 with 1, y\, y 2 and yij/2 
harmonic. Denote its transition probability kernel from time —m to time —n > —m by 
K~ m - n (y,dv). Then K~ m - n and K~ m - n are related via 

/(£) K~ m '- n (e, dx) = h{9) I ({o cf,- 1 ) (y) k- m '~ n (0(9),dy) V/ measurable. 

./[0,l)2V/i / 



'[0,00) 2 

In particular, 



[0,oo)2 

h{9) 



9i (x)K^°(9,dx) 



(F C N 9) 



(5.4) 

0< N<j, i = 1,2, 



'[o,i) 2 v ^ 

since {F 3 c g)i(0) = E[(F* g)i(X(-N))\X(-j) = 9} for all < A < j. For j G N, if we let 

f« = (y©(_ n )) neNo (5.5) 

and for all — n < —j 



denote the Markov chain Y started at time — j with Y^\—j) 
set yW(-n) = 0(0), then we can rewrite (15. 4p as 



(Fc9)i0) = H9)E [(§ o r 1 ) (^ 0) (0) 



h(9)E 



(F c N g)i 
h 



r A (yO)(-A)) 



(5.6) 



To establish (I2.25j) . and hence Theorem 12. 161 for c n = c, we need the following lemma, the 
proof of which is postponed. 

Lemma 5.1 For any fixed N G No, all weak limit points of {Y^\—N)}j^ as j — > oo are 
supported on (j){Roo) U ([0, 1) x {0}) U ({0} x [0, 1)). 
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We first complete the proof subject to Lemma 15.11 Without loss of generality, take i = 1. 
Note that, since g G TCq, we have g\{x) + g 2 {x) < C(l + xi)(l + £2) for some C > 0. 
Consequently, by the moment equations ()A.2[) ()A.3() . the family of functions 

(5.7) 



h 



fceN ,i?e[o,i) 2 



is uniformly bounded. Now fix 9 G [0,oo) 2 . If {i^}meN is any subsequence along which 
lim m _ +00 (Fc m 5)i(^) exists, then we can find a further subsequence {Jm}meN such that Y^ m > 
converges weakly to a limit Y°° = (Y°°(— n)) ng ^ as ([0, l] 2 ) N -valued random variables with 
the product topology. In particular, Y^ m \—N) converges weakly to N) for each iV G 

N . 

By Theorem 12.61 the family 



1 (?/) > (5- 



is continuous on [0, l) 2 . In fact, it is also continuous at (j)(Roo) with 
(F c k g) 



\2 



-o<frAffi = \ lt g VkeN , ze^Roc). (5.9) 



Indeed, this follows from these observations: (1) g G Hq, and hence ((<7i//i)°0 )(z) = \\z for 
z G 0(i?oo) and is continuous at z; (2) by (F* g)i{9) /h(9) = E[(gi/h) o ^-^(^(O))]; 

(3) because 3^ , i = 1,2, are martingales while 4>(Roo) = {(1,0), (0, 1), (1, 1)} are extremal 
in [0, l] 2 , it follows from the Markov inequality that K~ k,0 ((p(9),dy) converges weakly to the 
point mass at z as <fi{9) — ► z for z G (p^R^). By Lemma 15-H we can now substitute j m for j 
in (|5.6p and take the limit m — > oo, to obtain 



lim (F^g)^) = h{9) E 



ViVGN . (5.10) 



Denote the distribution of Y°°(— N) by fijsr. Again by Lemma 15. 1\ [in is concentrated on 
0(-Roo) U [0, 1) x {0} U {0} x [0, 1). Consequently, because ((F^g) 1 /h o vanishes on 

{0} x [0, 1], we have 

lim {Ftg)i{9) 

m— >oo 

= h{9) {^ N {{1, 1)}A 1)(00j00) + jf 1 o (yi, 0) m (d yi x {0})) . 

Since yi, 7/2, 2/12/2 are bounded continuous functions on [0, l] 2 and since 

and E[Y/ im) (-iV)y 2 (im) (-iV)] = <k(&)<h(.&) with <A = (<h.,<h), we mus t also have / yifi N {dy) 
= 4>i(9) and J yiy2l^N(dy) = <j)i(&)4> 2 {9). By our property of the support of /x/v> we thus find 

Miv{(l,l)} = ^(9)^0) = (5-12) 



ViHN(dyix{0}) = I yi(l-y 2 )m(dy)= 0i(0)(l - 2 (#)) = 77V- (5-13) 

/i(0) 
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Therefore 



lim {F 3 c m 9)1(6) - A 1)(oO)Oo) 0i0 2 - Ai i(oOiO) 0i 

< sup (^lio0- 1 )( 2/1 ,O)-A li(oOjO) y 1 

1/1 G [0,1] 

{Ffg^x, 0) - A li(oo>0) x 



n(6>) sup 

a;>0 



1 + X 



(5.14) 



Next, note that {{F^ g)\)N<=$ restricted to (0, 00) x {0} are the iterates of the renormalization 
transformation acting on diffusion functions on the halfline with initial diffusion function 
gi(x,0). Since ]im x — too gi(x,Q)/x = Ai 0) £ [0, 00), Theorem 5 of Baillon, Clement, Greven 
and den Hollander [3] implies that sup^Q \(F^ g)i (x, 0) — A 1 ( OOi0 )x|/(l + x) — > as N — > 00. 
(The case Ai^oo^) — is not included in Theorem 5 in [3], but an examination of the proof 
shows that the same result holds.) Since N can be taken arbitrarily large in (|5.14j) . we have 
established the convergence in (|2.25p along the subsequence {j m }meN- Since {(F£ g)i(9)}j & ^ 
is uniformly bounded, (12.250 now follows and the proof of Theorem 1 2 . 1 61 for c n = c is complete. 



We now prove Lemma l5.11 

Proof of Lemma 15.11 We must prove that the weak limit of {Y^ m ^} m ^, written Y°°, 
satisfies 

P(y°°(-JV) G (p{Roo) U [0,1) x {0}U{0} x [0,1)) = 1 ViVGN . (5.15) 
The proof consists of the following three steps: 
(A) Show that (Y°°(— n)) nS N ; * = 1, 2, are backward martingales on [0, 1], i.e., 



E 



YTi-k) <yr(-n)) 



n>k+l 



Y°°(-k-l) 



. = 1,2, 



(5.16) 



implying that lim^^oo Y°°(— n) = Y°°(— oo) exists a.s. by the backward martingale 
convergence theorem (see e.g. Section 4.6 in Durrett |18j). 



(B) Show that P{Y°°(-cx)) G 0(Roc) U [0, 1) x {0} U {0} x [0, 1)} = 1. 

(C) Show that P{y°°(-A^) £ (/)(Roo) U [0, 1) x {0} U {0} x [0, 1)} = 1 for all N G N . 

Since (Y^ m \— n)) ng N , m G N, i = 1, 2, are bounded backward martingale sequences, (A) 
follows from a general result on weak limits of backward martingale sequences, which we state 
as Lemma 15.21 below. The proof of (B) given below uses Lemma 14.81 which relies on uniform 
lower and upper bounds on {F™g} ne ^ , where assumptions (12.240 and g G TLq are crucial. 
The proof of (C) given below is achieved after approximating Y°° by the Markov chains Y^ m ' 

and using the fact that Y i , i = 1,2, are martingales. Note that it is not clear if Y°° is 
a Markov chain, because Y^ m ^ take values in [0, l) 2 while Y°° takes values in [0, l] 2 . Even 
though the transition kernels of y( Jm ) are consistent for m sufficiently large, they may not be 
(weakly) continuously extendable to [0, 1] 2 \[0, l) 2 . 
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Lemma 5.2 [Weak limits of backward martingales] 

For j G N, let Z^> = (Z^ (—n)) n ^ be a backward martingale, i.e., 



E 



zV\-k) {Z^\-n)) n 



>fe+i 



Z^X-k-l). 



(5.17) 



If {Z^\0)}j e fq are uniformly integrable, and Z^> converges weakly to a random variable 
Z°° = (Z°°(— n)) ng N in the space K N with the product topology, then n)) ng N is also 

a backward martingale. 



Proof. Since {Z^(0)}j^ o are uniformly integrable, we have 



Ve > 0, 3 N > such that E 
which is easily seen to be equivalent to 

Ve > 0, 3 N > such that E 



l^ (i) (0)|l|zu)(0)| 



>N 



< e 



[\ZW(0)\- N y 



< e 



V j G N, 



yj en. 



(5.18) 



(5.19) 



Since f{x) = (\x\ — N) + is a convex function, for all j, k £ N we have, by Jensen's inequality, 



E 



(\Z^(-k)\-Ny 



= E 

= E 

< E 

= E 



f(z®(-kj) 

f(E[zW(0)\(zW(-n))n>k) 



E[/(Z^(0))|(Z^(-n)) n > fc ] 
/(Z&')(0))1 =e\(\Z®(0)\-N) 



Therefore {Z^\— ^)}jeN,neN i s a uniformly integrable family. 

For each k E No and j G N, and any bounded continuous function /: 
martingale property of Z^> implies that 



E 



/ ((^ } (--))„> fc+ i) (z^{-k) - Z®(-k - 1)) 



0. 



(5.20) 



I, the 



(5.21) 



since converges weakly to Z°°, and {Z 3 {— k)}j^\ and {Z 3 (—k — 1)} jG n are uniformly 
integrable, we may pass to the limit j — > oo and obtain 



E 



>k+l 



Z°°{-k)- Z^i-k-l) 



0. 



(5.22) 



Indeed, the latter is easily verified by applying Skorohod's representation theorem, which 
allows for a coupling between {Z^}j e ^ and Z°° such that the convergence is a.s. From (I5.22j) 
we have 



E 



>k+l 



E 



Z°°(-k)-Z°°(-k-l) {Z°°{-n)) n>k+l 



0. 



which implies that 



E[Z°°{-k) - Z c 



-k-1) Z c 



- n ))n>k+l] = a - S - 



(5.23) 



(5.24) 



and thus establishes the martingale property for Z°°. 
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We are now ready to verify (B) and (C). 

(B) : Note that 

(K^) U ([0, 1) x {0}) U ({0} x [0, 1)) = ([0, 1] x {0}) U ({0} x [0, 1]) U (1, 1). (5.25) 

Suppose that (B) fails. Then there exists a u G (0, 1] \(1, 1) in the support of the distribution 
of Y°°(— oo). In particular, for each e > there exist 5(e) > and N(e) > such that 

p{y°°(-n) G B e/2 (u) Vra > N(e)\ > 5(e), (5.26) 

where B e / 2 (u) = {y G [0, l] 2 : \\y — u\\ < e/2}. Since converges weakly to Y°° as m —* oo, 

for each M G N we can find an m* = m*(M) sufficiently large such that 

p jf(im*)(_ n ) G Be (u) n (0,1) 2 ViV(e) < n < N(e) + Af} > -5(e). (5.27) 

We now derive a contradiction with Lemma [4.81 as follows. By assumption (|2.24p and the fact 
that g G 7i r , implying g\(x) + g 2 (x) < C(X + xi)(l + x 2 ) for some < C = C(g) < oo, F'^g 
satisfy the same upper and lower bounds for all n G N. It is then easy to check that in Lemma 
IMlwith A = {F™g} ne n Q condition (|4.17p is satisfied for all G (0, oo) 2 , and conditions (|4.19l - 
I4.2UP are satisfied for all a > (and the analogue of (|4.19|) for vertical strips). Since the 
transition kernel K~ n ~ ~ n ((f)(9), dy) is related to the biased equilibrium measure TJ c (dx) 

through the coordinate change (j), Lemma I4.8I (i) and (ii) imply that, for u G (0, 1] 2 \(1, 1) and 
e > sufficiently small, 

mf ^"^-"(^[O,!) 2 ^^)) > 0. (5.28) 

S6B f (S)n[l),l) 2 

This uniform rate of escape from B e (u) contradicts (|5.27p . where M can be chosen to be 
arbitrarily large while 5(e) > remains fixed. 

(C) : For e > 0, let 

U £ = \y G [0,1] 2 : inf \\y - z\\ < e\ . (5.29) 

I ze<t>(Roo )u[o,i)x{o}u{o}x [o,i) J 

Since Ymin^^Y 00 (—n) = Y°°(— oo) a.s., we can choose M = M(e) sufficiently large such 
that F(Y°°(-M) G U e ) > 1 - e. Since Y^ m \-M) -> Y"°°(-M) in distribution as m -> oo, 
we can choose m*(M) sufficiently large such that F(Y^ m \—M) G U 2e ) > 1 — 2e for all 
m > m*. By the geometry of <j)(Roo) U ([0, 1) x {0}) U ({0} x [0, 1)) and the fact that Y} jm) (-), 
i = 1,2, are martingales for the Markov chain (Y^ m \— n)) n6 N , an elementary application of 
the Chebychev inequality shows that, for all m > m* and L > 2, 

P (yO™) (0) G U 2Lt ) > (1 - 2e) (l - |) . (5.30) 

By the weak convergence of y^ m ^(0) to y°°(0) as m — > oo, the same holds for Y°°(0). Now 
let e — > and L — > oo such that eL — > 0. Then we find that 

P (y°°(0) G 0(«oo) U ([0, 1) x {0}) U ({0} x [0, 1))) = 1. 
The same argument works for Y°°(—N) for any iV G No- ■ 
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6 Proof of Theorem 12.161 with varying c n 

Proof. The proof of Theorem 12.161 with varying c n follows the same line of argument as that 
for constant c n , except for a few technical differences, which we now outline. For the rest of 
the section, let (X (—n)) n ^ denote the backward time-inhomogeneous Markov chain with 
transition kernels 

¥(X(-n) G dx I X(-n - 1) = 9) = T ?^ 9 {dx), (6.1) 

and let {X h (—n)) n ^ denote X /i-transformed by h(x) = (1 + x\)(l + x 2 ), which is still a 
harmonic function for X. Both X and X h generalize their counterparts in Section [5j We 
proceed by first establishing the analogue of Lemma 15.11 where {^-^bjeN are now defined in 
terms of our current X and X h . 

The proof of Lemma 15.11 in Section [5] is based on Lemma 14.81 which no longer applies in 
our current context, because if c n can be arbitrarily large, then we lose the uniformity of the 
escape probability with respect to {I^' } nS N - So, the first task is to formulate a suitable 
analogue of Lemma [4.81 for our current X and X h , which would imply the analogue of Lemma 
15.11 for the present context. In the derivation of Theorem 12.161 for constant c n from Lemma 
15. 1| we used the following fact from Baillon, Clement, Greven and den Hollander [HJ: for the 
renormalization transformation F c acting on one-dimensional diffusion functions / : [0, oo) — > 
[0, oo), where / is positive and continuous on (0, oo), locally Lipschitz at 0, /(0) = and 
lirn^oo f(x)/x = A G [0, oo), we have sup x>0 \{F™f){x) — As|/(1 + x) — > as n — > oo. Our 
second task is therefore to establish the analogous result for F^f. The two technical points 
outlined above will be addressed in Lemma 16. II and Proposition 16.21 below. 

Observe that, by Proposition IA.lt f° r au ~ m < ~ n 5; an d 9 G [0, oo) 2 , the backward 
Markov chain X satisfies the moment equations 

E[X(-n)\X(-m) = 9\ = 9, (6.2) 

E[X 1 (-n)X 2 (-n) | X(-m) = 9\ = 9 X 9 2 , (6.3) 

E[X,(-n) 2 | X(-m) = 9\ = 0? + I V) - I (F^g) t (9), i = 1,2. (6.4) 




E[ 9i (X(0)) | X(-m) = 9\= (F^gU9), i = 1,2, (6.5) 

From the point of view of variance increment, (j6.4j) indicates that the natural time associated 
with (X(— n)) n >o is not n, but rather Y17=o • Therefore to obtain a uniform bound on 
escape probabilities for the Markov chain X h , we formulate the analogue of Lemma 14.81 as 
follows. 

Lemma 6.1 [Uniform rate of escape of {X h {— n)) n >o from small balls and thin strips] 

Let (c n ) nG N an d 9 be as in Theorem \2.1b\ Let (X(—n)) ne fq denote the inhomogeneous back- 
ward Markov chain with transition kernel (jg.jp . and let (X h (—n)) n ^jq denote (X(—n)) n ^ 
h-transformed by h(x) = (1 + xi)(l + x 2 ). There exists an increasing sequence (nk)k£N C No 
with n = such that J]^* 1 " c i l G I^" 1 ' M f or some A > 1 /or all k G N . For A C [0, oo) 2 , 
denote r7 m = inf{-j > -m : X h {-j) £ A}. Then 



(i) For each 9 G (0, oo) 2 , there exists e > such that 

m f tor * -** I = *) > o- ( 6 - 6 ) 

sgs e (e) 
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(m) For each a > 0, there exist e, N > such that 



inf 1 

fe£N 

x£E[iV,oo) X [ck — e,a + e] 

inf 

fcen 

cc6[a — e,c*+e]x [iV,oo) 



T ~ nk+1 < -nu 

1 [N,oo)x[a-e,a+e] — nk 



_-"fc+i 

[a— e,a+e] X [N,co) 



X h {-n k+l ) = xj > 0, 
* h (-n*+i) = f) > 0. 



(6.7) 
(6.8) 



Proof. The existence of the increasing sequence (n&)fc e N with the prescribed property follows 
immediately from our assumptions that inf ng N c n > and ^ngN c n ] 



oo. The rest of the 



proof parallels that of Lemma 14.81 First we prove (jO|) - (pT8|) with X h replaced by X. By 



(|4.22p . for each m G N and f € [0, oo) 



1,2, conditioned on X(— m — 1) = x, we have 



E 



L(l + X 4 (-m))2 



1 



> 



1 + Xi 

1 



-E 



1 



l + Xi(-m) 
2 



+ 



E 



(FM g ).(x(-m)) 
(l + X i (-m))3 



(1+Xi) 2 C m (l+Xi) 



E 



(f [ml g),(I(-m)) 
(l + Xi(-m)) 3 



(6.9) 



where we applied Jensen's inequality. Conditioned on X(— rik+i) = x, we can apply (|6,9 
iteratively to obtain, for i = 1, 2, 



E 



1 



(l + X(-n fc )) 2 



> 



(6.10) 



"fc+i— 1 



(1 +x. 



^2+ E 



-E 



(F [ml ff),(I(-m)) 



l + JSQ(-m-l) (l + Xi(-m)) 3 



If (16. 6p fails when X' 1 is replaced by X, then there exists 9 G (0, oo) 2 such that, for all e > 0, 
there exist sequences few j oo and a?w g 2? e (0) (depending on e) such that 



lim P[r 



t(0 H 



B e (0) 



1 < ~ n kW 



X(-n kW+l ) = x 



(I) 



0. 



(6.11) 



Now we apply (|6.10p to X{—n) for — n k (i) +1 < —n < —n k {i) with X(— n k (i) +1 ) = x*w. By 
(16. lip , as I — > oo, the two sides of (I6.10P satisfy 



l.h.s. < 



r.h.s. > 



1 



(1 + 0i - e 
1 

(1 + 0i + e) 2 



|2+o(l), 

+ E 



2 , , „ „ 5 
— l-o 1 

Cm (1 + 0t + e) 4 



(6.12) 
(6.13) 



where in ()6.13p we have used the assumption that gi(x) > a%x% + j3iX\X2 for some aj,/3j > 
and oti + fa > 0, i = 1, 2, which implies that {i ? ' m '<7}meNo satisfy the same lower bound and 
(fN g )j(f) > 5 > uniformly for x G £ e (0) and m G N . Since Em=n fc ^ > A" 1 > 
uniformly for all k G No, (|6.12p and (|6.13p are incompatible for e > sufficiently small and 
I G N sufficiently large. Therefore (|6.6[) must hold for X in place of X h . The proof of 
(|6.7p - (|6.8p for X in place of X is similar, and we leave the details to the reader. 

To verify that (|6.6p also holds for X h , we apply Lemma 14.51 and note that h{x) = (1 + 
X\){\+X2) is bounded uniformly from above for x G B e (6), and bounded uniformly from below 
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by 1 for x E [0, oo) . The proof of (|6.7p - (|6.8p for X is essentially the same as its counterpart 



in the proof of Lemma [4.81 Note that by Lemma [4.51 the law of X h (t< n ^ x , a _ e a+e ^ A (—rife)) 
conditioned on X h {— rik+i) = x € [N, oo) x [a — e, a + e] is absolutely continuous with respect 



to the law of X(t, 
h(-) 



-"fe+i 

[JV,oo)x [a— e,a+e] 



A (— rtfc)) conditioned on X(— n k+ i) = x, where the density 
is ^gy- As in the proof of Lemma 14.81 it suffices to show that for any fixed < e < a < oo, 



lim sup 

x 2 ^ [ a ~ e,Q!+e] 



X(-n k+ i) = x)=0. (6.14) 
Since (X±(— n)) n < nk+1 is a martingale, by Doob's inequality and (|6.2H6.4p . we have 



r [a5)xR - nk 



X(-n k+ i) = x 



< 



x\ 



sup \X\(-n) -x\\ > — X(-n k+1 ) 
-rit-Li <— n<— m, ^ 



< 5 = 



(Xi(-n fc ) - xi) 2 X(-n k+1 ) = x 



16 



1 \ n=n fc Cn J 



(6.15) 



Note that ^n=~n fc 1 5~ — A uniformly in fc. Since <? £ Hq, we have <?i(x) + 52 (aQ < ^(1 + 
x\)(l + X2) for some K E (0, 00), and by Proposition IA.ll {F^g} n ^ all share the same 
upper bound. Equation (|6.14|) then follows immediately. I 

Remark. Note that (|6.6p - (l6,8p with X in place of X h are proved using only the assumptions 
that X^neNo c n 1 = 00 an d' {-^ n 'fl , }neN have a uniform lower bound which is positive and 
uniformly bounded away from on (a, oo) 2 for each a > 0. Only in deriving (|6.7|) - (|6.8|) 
from their analogues for X, did we use the assumptions that inf ne jjo c n > and, {F^g} n£ ^ 
have a uniform upper bound (p = (<f>\, <f>2 ), where 4>\{xx,X2) grows sub-quadratically in x\ and 
(^2(^11^2) grows sub-quadratically in X2- 

Using Lemma 16.11 and the fact that 1, x%, X2, x\x% are still harmonic functions for the 
Markov chain {X(—n)} n ^ , we deduce the analogue of Lemma 15. II in our present context by 
the same arguments as in the original proof. To deduce Theorem 12.161 with varying c n from 
the analogue of Lemma 15.1} we need to address the second technical point outlined at the 
beginning of this section. 

Proposition 6.2 [Convergence to fixed points under fN : the half line case] 

Let (c n ) n£ N satisfy X] ngNo c n l = 00 • Let f(x) : [0, 00) — > [0, 00) be positive and continuous 
on (0,oo), locally Lipschitz at 0, /(0) = and lim^oo x~ 1 f(x) = A G [0, 00). Then we have 



lim sup 



x>0 



(F^f)(x)-\x 



1 + x 



0, 



(6.16) 



where F<- n > are renormalization transformations acting on one- dimensional diffusion functions. 

Proof. Note that we do not require inf ng N c n > as in Lemma 16.11 The case c n = c is 
covered by Theorem 5 of Baillon, Clement, Greven and den Hollander [3]. Here we give a 
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proof along the same line of argument as we have been pursuing so far in this section for the 
proof of Theorem 12.161 with varying c n , except that we do not need, to appeal to the current 
proposition. 

As in Section 2.3 of Baillon, Clement, Greven and den Hollander [3J, we make use of the 
concave upper envelope / + and the convex lower envelope /~ of /. It is easy to see that / + and, 
/~ in the case A > 0, satisfy the same constraints as specified for / in the proposition. Since, 
for any c > 0, F c is convexity preserving and order preserving by Proposition 3 of [3 J , together 
with Jensen's inequality we have, for each x G [0,oo), (F^f+){x) | f+(x) and j 
f~{x) for some /+ and /" as n -> oo, and (F^f~)(x) < (F^f)(x) < (F^f+)(x) for all 
n G No- We claim that it suffices to show that /^(x) = /^(x) = Ax. Indeed, 

(l-y) (F^f)o^-\y)-\y\, 
where ^i(x) = j^;- Since F c preserves the slope at infinity, we have that 

My) = (i-y)(F [n] f)°^Hy), (6.17) 

^-(«) = (l-y) (F^r)o^Hy), 

are all continuous functions on [0, 1]. If /^(x) = /^,(x) = Ax, then, on [0, 1], Vn iv) decreases 
monotonically to Ay as n —> oo, while ipniv) increases monotonically to Ay as n — > oo. 
Since the monotone convergence of a sequence of continuous functions to a continuous limit 
is necessarily uniform on compacts, the sup-norm convergence of ip n (y) to Xy on [0, 1] follows 
since ip n is sandwiched between and ip~ . 

The proof that /+(x) = lim„_ 00 (FN/+)( x ) = Ax and f~(x) = )hn^ 00 {FW f~)(x) = Ax 
now follows the same argument as that used for Theorem 12.161 with varying c n . First consider 
the case {i^ n V + }neN W1 th A > 0. In the proof of Theorem 12. 161 with varying c n , we replace 
X there by the [0, oo)-valued Markov chain {X{— n)) ng N with transition kernels 

P(X(-n) G • | X(-n - 1) = x) = r^ W / + (•); 

X h is replaced by X h , which is the /i-transform of X by the harmonic function h(x) = 1 + x; 
4>{x) is replaced by 4>\{x) = in Lemma l5TT| the relevant boundary points now consist of 
only {0} U {1}. Lastly, because of the one-dimensional setting, we only need to establish the 
analogue of (16. 6p . By the remark following the proof of Lemma 16.1} the only assumptions we 
need here are ^neN c n l = 00 anc ^' a uniform lower bound on {i ? f n '/ + } n eN which is positive 
and bounded away from on [a, oo) for each a > 0. Note that / _ provides such a lower 
bound. The case {F™ f~} ne ?$ with A > is identical. For the case A = 0, we only need to 
consider {i ? ^/ + }neN - Everything remains the same, except that the uniform lower bound 
on {F[ n ]/ + } ng N is now provided by /+. Indeed, as a limit of concave functions, /+ is also 
concave, hence either /+ = 0, in which case we are done, or /+ is positive and non-decreasing 
on (0, oo), which is sufficient for the proof of the analogue of (|6.6p to go through. I 

With Lemma [6. II and Proposition [621 we can now proceed as in the proof of Theorem 12. 161 
for constant c n and extend it to varying c n . We leave the details to the reader. I 



sup 

x>0 



- Ax 



1 + x 



sup 

2/6(0,1) 
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A Appendix 1: Moment equations and estimates 



Proposition A.l [Moment equations and estimates] 

Letg eC,6 £ [0, oo) 2 , c > and let T C J 5 be any equilibrium distribution of (\2.1\i with generator 
0Q|. Let X = (Xi,X 2 ) be a random variable with distribution T^ 9 . Then: 

u 

(i) For any f(x) E C? ([0, oo) 2 ) that differs from a function with compact support by only a 
constant, 



E' 



'-■•9 



(L c /f)(X) 



E 



1=1 ' 1=1 1 



0. (A.l) 



(ii) For all g G TC a with < a < c, all 9 £ [0, oo) 2 and i = 1, 2, 

E% a [Xi] = 0i, (A.2) 
E < 2 9 [X 1 X 2 ] = 9 1 9 2 , (A.3) 

Ef[Xf] = 9? + -E^X)] = Q ^ + \^9W), (A.4) 
where all expectations are finite, 
{in) Let g € TL a with < a < c, and let K be any compact subset of [0, oo) 2 . Then 

sup Ef 9 [(Xi + X 2 + 2) 2 log(Xi + X 2 + 2)] < (A.5) 

d>cfi&K 

for some C Ct x,g < °o depending only on c, K and g. Consequently, g\ and g 2 are 
uniformly integrable with respect to {^' 9 } c iy c q^k~ 

Proof, (i) This part follows from the observation that, with our choice of /, 

f(X(t)) - /(X(0)) - J\l c / f)(X(s))ds (A.6) 

is a martingale. Taking expectation and noting the stationarity of the distribution of X(t), 
we obtain (|A.ip . 

(ii) We first prove that the expectations in (|A.2HA.4~|) are all finite. Once this is settled, the 
equalities will follow easily. 

Finiteness : Let h 6 C 2 ([0,oo)) be such that h(r) = r for r G [0,1], h is constant on [3, oo), 
ti G [0,1] and h" G [-1,0]. Let h n {r) = nh(^). Then h' n G [0,1], K G [~^,% and h n {r) | r, 
h' n (r) | 1, ^(r) — ► as n — > oo. 

(|A.2j) : We apply (1A.1|) for f(x±,x 2 ) = h n (pix\ + ^2^2) with fixed pi,p 2 > 0. Since (in the 
formulas below we suppress the argument) 

d Xl h n (p lXl + p 2 x 2 ) = pih' n , ^ 
d 2 h n (pixi + p 2 x 2 ) = p 2 h" n , 
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and h n {p\X\ + P2X2) differs from a function with compact support by a constant, by substi- 
tuting the partials into (|A.1|) . we get 



c,9 



E 



which can be rewritten as 



i=l 



cE 



(piXi + p 2 X 2 )h' n 



= cE c 4 9 
9 

< cE c J 9 



(pi0i + P20 2 )ti n 

Pi 01 + P202 



+ E' 



'-■■9 



0, 



(PlSl + 0202 



(Ai 



(A.9) 



since /i" < and 51,52 > 0. By monotone convergence as n — > 00, we get 

PiE^Xi] + p 2 E^[X 2 ] < p^x + p 2 fl 2 . 

Since pi,p 2 are arbitrary, we obtain E^.' 9 [Aj] < 0« < 00, i = 1,2. 



(A.10) 



()A.3[) : Here we apply (jA.ip for /(xi, x 2 ) = /i n ((l + £i)(l + x 2 )) . The calculations are similar 
to that for (|A.2|) . which we skip. 



(|A.4|) : Here we apply (jA.ip for f{x\,x 2 ) = h n (pix1 + p 2 x|) with fixed pi,p 2 > 0. Since 



d Xi h n (pix\ + p 2 x\) = 2piXih' n , 
d 2 Xi h n {pix\ + p 2 x\) = 2 Pi h' n + ipfxjh'n, 

by substituting the partials into (jA.ip . we get 
E^[2cp 1 Xi(0 1 -X 1 )/4 + 2cp 2 X 2 (0 2 -^ 
Rearranging terms, we obtain 



(A.ll) 



0. 



2c E' 



eg 



( Pl xf + P2 xi)h' n 

2cE c J s 
9 

+ E C J 9 



(piOyXx + p 2 02^ 2 )^n] + [(Pl5l + P29 2 )h' n 

\{p\Xlg l + p\xlg 2 )U' n 
< 2c { Pl 6l + p 2 9l) + 2E<2 9 Up l9l + p 2 g 2 )h' n 



(A.12) 



Since g S TL a with < a < c, we have g\{x) + g 2 {x) < C(l + a?i)(l + x 2 ) + a(x\ + x%). 

Substituting this bound into (|A.12[) and setting p\ = p 2 = 1, using the fact that E^ 9 pQ] < 6i 

9 

and [XiX 2 ] < 00, and rearranging terms, we get 



< C < 00. 



(A.13) 



2{c -a)E c f[{Xl + Xl)h' n 
By monotone convergence as n — > 00, we obtain El' 9 ^?] < 00. This also implies E 6 .' 9 ^] < 00. 

9 9 

Equality: Having thus proved that the expectations in (|A.2I IA.4|) are finite, we are now ready 
to prove that equality holds. To that end, return to (|A.9p . Since E C g 9 [p\g\ + p|<72] < 00, fe" € 
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[— —j 0] and h'n —* as n — > oo, (|A.2p follows by applying the dominated convergence theorem. 
By the same argument, (|A.4|) follows by applying the dominated convergence theorem to 
(|ATl2|> . provided that 



(/ofxfgi + P^252)K(pix? + ,02x1)1 < C(x? + x|) (A.14) 

for some C < oo independent of n. To see the latter, note that h'^(r) = ^h"(^) = -h"(-)- < 
since h" € [-1,0] and h"(^) ^ only when £ < 3. The bound in fOIj) then follows 
readily. 

To verify (|A.3j) . we apply (lA.lj) for f{x\,x 2 ) = ^((^l+a^) 2 ) instead of /i n ((l+xi)(l+x 2 )). 
This gives 



E' 



'-■9 



2c(0i + 9 2 - X x - X 2 )(X 1 +X 2 )h' n + 2( 5 i + 52 )< + A(X 1 + X 2 ) 2 ( 5l + g 2 )h!' n =0. (A.15) 



Since (xi + x 2 ) 2 ^n(( x i + x 2 ) 2 ) < 3 and ¥, c 4 g [gi + g 2 ] < 00, we can apply the dominated 
convergence theorem in (|A.15P as n — > 00. Then, together with (|A,2p and (|A.4p . we obtain 
(POD . 

(hi) This part follows from similar computations as in part (ii). Let c' > c be arbitrary, and 
abbreviate Xi = 1 + Xi, 0{ = 1 + 9i, x% = 1 + X{ for i = 1, 2. We first show that 

E d J 9 \X x X 2 log(X! + X 2 )] < 00 (A.16) 

by applying (jA.lh to h n {x\X2 \og{x\ +x 2 )) . Then we apply (jA.lh to /i n ((xi +x 2 ) 2 log(Si +x 2 )) 
to prove (|A.5|) . 

(|A.16h : Let f(xi,x 2 ) = h n (xix 2 log(xi + x 2 )), which differs from a function with compact 
support by a constant. Since 



x 2 log(xi + x 2 ) + r — — — h' n , 



Xl+X 2 . 

x 2 x\ \ , /_ XiX 2 \ 2 , 



3 /t n (xix 2 log(:ci + x 2 )) = (- — — - + -n — 2 .„ 1/4 + (x 2 log(xi + x 2 ) + - — — ^) /4, 

and since the same holds if we interchange the indices 1 and 2, by substituting the partials 
into (jA.ip and noting that d^hn < 2h' n ,d^ 2 h n < 2h' n , we get 



E** c'fr - Xi)(x 2 log(X 1 + X 2 ) + 3^L- ) h' n 

+ c'(9 2 - X 2 ) (X X log(Xi + X 2 ) + * lX \ ) h' n + 2( 51 + g 2 )h' n ~ 
V Ax + A2 / 

Rearranging terms and noting that %^+j£ 2 < X\ A X 2 , vre find that 
2c'Ef 9 [X 1 X 2 log(X 1 + X 2 )ti n ] < Ef 9 [c , (9 1 X 2 + X 1 9 2 )(l+log(X 1 +X 2 )) + 2( 9l +g 2 )\. (A.18) 

By assumption, g\{x) + g 2 (x) < C(l + xi)(l + x 2 ) + a(x\ + x 2 ). Substituting this bound 
into (IA.18p . applying monotone convergence as n — > 00, and noting that (IA.4P implies that 
E 6 .' 9 ^? + A"|] < (f>(0i,6 2 ) for some quadratic polynomial depending only on c and 5, we 
easily verify that 

Ej 9 [JM 2 log(Xi + X 2 )] < 0(0!, 2 ) (A.19) 



> 0. 
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for some cubic polynomial depending only on c and g. 

By applying (jA.ip to h n ((x 1 + X2) 2 log(xi + £2)) and using (|A.19p . it can be shown that 

^f[{X x + X 2 f log(*i + X 2 )] < ^1, 5 2 ) (A.20) 

for some cubic polynomial <fi depending only on c and g. The uniform bound in (IA.5P then 
follows. The calculations, which we omit, are similar as before. 

Since gi(x) + g 2 (x) < C{x\ + x 2 ) for some C < 00, which by (|A.20|) is uniformly integrable 
with respect to {r^' 9 } c / >c g eK for any compact K C [0, oo) 2 , it follows that g\ and g 2 are also 
uniformly integrable. I 

Remark. By similar computations, it can be shown that (IA.5P is still valid when the logarithm 
in the left-hand side of the inequality is raised to an arbitrary power. 



B Appendix 2: Properties of uniformly elliptic diffusions 

In this Appendix, we list some facts about uniformly elliptic diffusions that are needed in the 
proof of Theorem 12. 31 We thank S.R.S. Varadhan for pointing out some of the relevant results 
and references on uniformly elliptic diffusions. 

Theorem B.l [Uniformly elliptic diffusions in M. d ] 

Let b : M. d — > M. d be a bounded measurable map, and let a : M. d — > be a continuous map, 
where Sd is the space of symmetric non-negative definite d x d real matrices. Assume further 
that a(-) is uniformly elliptic, i.e., there exists < A < 00 such that for all f,fl € M. d , 6^0, 

A -i<<Mj^ <A , 



Then, for each x € M. d , the martingale problem with generat 



or 



Lf=J2 W&g^r.fW +X>(*)^:/(2), / 6 C!(R d ), (B.l) 

i,j=l 1 ^ i=l 1 



has a unique solution W x in the space of probability measures on £1 = C([0, 00), ~R ) with 
F x (lu G f2 : lo(0) = x) = 1. The family of solutions {^ x }xeR d defines a strong Feller and strong 
Markov process that admits a transition probability density pt(x,y) with respect to Lebesgue 
measure for each t > and x G M. d . Furthermore, for each t > and x* € M. d , 

lim \\p t {x, •) -p t (x*, -)||i = Jim / \p t {x, y) - p t {x* , y)\dy = 0. 

Proof. All facts follow from results in Stroock and Varadhan [33J. For the well-posedness 
of the martingale problem, see Theorem 7.2.1 therein. For the strong Markov property, see 
Theorem 6.2.2. For the strong Feller property, see Theorem 7.2.4. For the existence of the 
transition density, see Theorem 9.1.9 and Lemma 9.2.2. Lastly, for the Li -continuity of the 
transition density, see Theorem 11.4.3. I 
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Theorem B.2 [Diffusions restricted to bounded domains] 

Let a and b satisfy the conditions in Theorem IB. 11 and let {F x }^ e Rd denote the family of 
solutions to the martingale problem with coefficients (a,b) in (IB.lj) . If a : ~R d —> Sd and 
b : M. d — > M. d are locally bounded measurable maps with a = a and b = b on a bounded open 
set D, then for any x G D and any solution V x to the martingale problem with coefficients 
(a,b), ¥ x = ¥ x on T TD , the sigma-field on Q generated by the family of projection maps 
{tt s : — » R d | it s (lu) = u(s A t d )} s >o> where t d (uj) = ird{t > : uo{t) £ D}. 

Proof. See Theorem 10.1.1 in Stroock and Varadhan 1 33 1 . 



Corollary B.3 [Transition density for diffusions restricted to bounded domains] 

Let b : M. d — > M. d be a locally bounded measurable map, and let a : W 1 — > Sd be continuous 
such that the martingale problem with coefficients a and b in (jB.ip is well-posed. Assume 
further that a is non- degenerate on D for a simply connected bounded open set D C M. d 
with smooth boundary. For any x G D, if F x is the solution of the martingale problem 
starting from x, then, for each t > 0, the measure fj,P(x,-) on Borel-measurable sets defined 
by n®(x,-) = F x (u> : t < T£)(uj),u>(t) G •) admits a density pP(x,y) with respect to Lebesgue 
measure. Furthermore, for each x* G D, there exist e, 5 > sufficiently small such that, for 
all x, x' G B e (x*), the ball of radius e centered at x* , the overlap between fJ,f(x, •) and f-tf(x', •) 
satisfies 

rf(x,D)+rf(x',D) - \\pf(x,-) -pf (gV)j|i 1 

2 - 2' [ ' ' 

Proof. By our assumptions on a, b and D, we can find coefficients (a, b) on R rf such that 
(a, b) = (a, b) on D, (a, b) are bounded, a is continuous and uniformly elliptic on R rf . For 
instance, we can define b = b on D and b = on M. d \D, define a = a on D and a = I on 
M. d \B where B is a large open ball containing D, and on B\D define d to be the harmonic 
interpolation between its values on dB and dD. By Theorem IB. H the martingale problem 
with coefficients (a,b) has a unique family of solutions {P x }^ eK d, which is strong Markov 
and admits a transition density pt(x,y) for all t > and x G M. d . By Theorem IB. 2\ for 
x G D, f x = F x on T TD . In particular, fJ,f(x,-) = ^(x,-) = f x {u : t < T D (u)),u(t) G •}. 
Since /^(x, •) is absolutely continuous with respect to ¥ x (u> : u(t) G •) with density pt(x,y), 
•) = Jif{x,-) also admits a density p^{x,y) with respect to Lebesgue measure for all 
x G D and t > 0. 

It is not difficult to see that the left-hand side of (|B.2|) is the mass of the maximal positive 
measure that is dominated by both fj,?(x,-) and fj,?(x',-). To verify (|B.2p . fix x* G D and 
choose e' > such that B2 e '(x*) C D. Then we can choose 5 > sufficiently small such that, 
for all x G B e i(x*), ¥ x (td < 5) < ^. To verify this claim, note that, given z G B € /(x*), if we 
define f(x) = \\x — z\\ 2 = Yli=i( x i ~ z i) 2 '■> then 

rtATD 

f(X(t A TD )) - f(X(0)) - / Lf{X{s))ds 

Jo 

is a martingale, where (X(s)) s >o has law F z . In particular, 

(e') 2 P z >D < S) < E[\\X(5At d )-z\\ 2 } (B.3) 



E 



/ 2^ ^(X( S ))(^(s) - *) + a«(A:(s))J 

^° i=l 



< <^ Cf),a,b, 
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where Co,a,b depends only on D and (a, b) on D. Therefore F z (td < 5) < 5CD, a ,b(^') 2 
uniformly for all z G B e i(x*). Choosing 5 sufficiently small, we then verify the claim. 

Applying Theorem IB. II to {P x }g^d, we can choose e G (0, e') small such that, for all x G 
B € (x*), \\ps(x, •) -ps(x*, -)lli < Jo, and hence, for all x,x' G B e (x*), \\pg(x, •) -ps(x', -)||i < g- 
Since for z G B e (x*), \\p s (z, -)-pf (z, -)||i = F(td < <5) < ±, we have (f, (£', -)||i < f 

for allx,x' G B e (x*). Finally, note that, for s,f G B e {x*), fJ<g(x,D) = 1-P s (t d <8)>l-\ 
and the same holds for fif(x',D), hence, substitution of all the estimates into the left-hand 
side of (|B.2I) yields the desired result. I 

Remark. Note that the constant on the right-hand side of (|B.2j) can be made arbitrarily 
close to 1 by choosing e, 5 sufficiently small. 

Theorem B.4 [Support theorem for uniformly elliptic diffusions] 

Let a,b,D and {P x }^£> be as in Corollary IB. 31 For any x G D, e > 0, and any continuous 
function ip : [0, t] — ► D with ip(0) = x, 

F s (u : sup \u(s) - ip(s)\ < e) > 0. 

Proof. The support theorem is a classic result of Stroock and Varadhan. The statement 
above follows Theorem (2.5) in Chapter V of Bass [4] and Theorem IB. 2 1 above. 

Theorem B.5 [Occupation time measure for uniformly elliptic diffusions] 

Let a, b, D and {P^jxG-D be as in Corollary IB. 31 If A C D has positive Lesbegue measure, 
then, for all x G D, K x [ fj 75 l^^g^GLs] > 0, where E, x denotes expectation with respect to ¥ x , 
and t d = inf{t > : u(t) (£ D}. 

Proof. The statement above follows from Theorem (8.5) in Chapter V of Bass [1] (which goes 
back to Krylov) in combination with the support theorem, Theorem IB.4| and the Girsanov 
transformation (see Theorem 7.2.2 in Stroock and Varadhan [33]). I 

Acknowledgment: The work in this paper was supported by DFG and NWO, as part of the 
Dutch-German Bilateral Research Group on "Mathematics of Random Spatial Models from 
Physics and Biology". AG and JMS were supported by the DFG-grant GR 876/12-1 - 12-3. 
DD was hosted by EURANDOM during two visits and is supported by an NSERC Discovery 
Grant. JMS is supported by GACR grant 201/06/1323. JMS and RS received travel support 
from the ESF scientific program "Random Dynamics in Spatially Extended Models" . FdH and 
RS are grateful to the Pacific Institute for the Mathematical Sciences and the Mathematics 
Department of the University of British Columbia, Vancouver, Canada, for hospitality: FdH 
from January to August 2006, RS from mid- April to mid-May 2006 when part of the work in 
this paper was completed. RS was a postdoc at EURANDOM from October 2004 to October 
2006. DD, FdH and RS thank Ed Perkins for valuable discussions. The authors thank the 
associate editor and the referee for an exceptionally careful reading of the paper and many 
helpful suggestions. 



45 



References 



[I] S.R. Athreya, M.T. Barlow, R.F. Bass and E.A. Perkins, Degenerate stochastic differential 
equations and super-Markov chains, Probab. Theory Relat. Fields 123 (2002) 484-520. 

[2] J.-B. Baillon, Ph. Clement, A. Greven and F. den Hollander, On the attracting orbit 
of a non-linear transformation arising from renormalization of hierarchically interacting 
diffusions, Part I: The compact case, Can. J. Math. 47 (1995) 3-27. 

[3] J.-B. Baillon, Ph. Clement, A. Greven and F. den Hollander, On the attracting orbit 
of a non-linear transformation arising from renormalization of hierarchically interacting 
diffusions, Part II: The non-compact case, J. Funct. Anal. 146 (1997) 236-298. 

[4] R.F. Bass, Diffusions and Elliptic Operators, Springer, New York, 1998. 

[5] R.F. Bass and E.A. Perkins, Countable systems of degenerate stochastic differential equa- 
tions with applications to super-Markov chains, Electron. J. Probab. 9 (2004) 634-673. 

[6] R.F. Bass and E.A. Perkins, Generalized mutually catalytic SDE's: uniqueness, in prepa- 
ration. 

[7] J.T. Cox, D.A. Dawson and A. Greven, Mutually Catalytic Super Branching Random 
Walks: Large Finite Systems and Renormalization Analysis, Memoirs of the American 
Mathematical Society 809, American Mathematical Society, Providence, Rhode Island, 
2004. 

[8] J.T. Cox, A. Greven, Ergodic theorems for infinite systems of locally interacting diffusions, 
Ann. Probab. 22 (1994) 833-853. 

[9] D.A. Dawson, L.G. Gorostiza and A. Wakolbinger, Degrees of transience and recurrence 
and hierarchical random walks, Potential Analysis 22 (2005) 305-350. 

[10] D.A. Dawson and A. Greven, Multiple scale analysis of interacting diffusions, Probab. 
Theory Relat. Fields 95 (1993) 467-508. 

[II] D.A. Dawson and A. Greven, Hierarchical models of interacting diffusions: Multiple time 
scales, phase transitions and cluster formation, Probab. Theory Relat. Fields 96 (1993) 
435-473. 

[12] D.A. Dawson and A. Greven, Multiple space-time analysis for interacting branching mod- 
els, Electron. J. Probab. 1 (1996), Paper no. 14, pp. 1-84. 

[13] D.A. Dawson, A. Greven and J. Vaillancourt, Equilibria and quasi-equilibria for infinite 
collections of interacting Fleming- Viot processes, Trans. Amer. Math. Soc. 347 (1995) 
2277-2360. 

[14] D.A. Dawson, A. Greven and I. Zahle, Continuum limits of multitype population models 
and renormalization, in preparation. 

[15] D.A. Dawson and P. March, Resolvent estimates for Fleming- Viot operators and unique- 
ness of solutions to related martingale problems, J. Funct. Anal. 132 (1995) 417-472. 

[16] D.A. Dawson and E.A. Perkins, Long-time behavior and coexistence in a mutually cat- 
alytic branching model, Ann. Probab. 26 (1998) 1088-1138. 



46 



[17] D.A. Dawson and E.A. Perkins, On the uniqueness problem for catalytic branching net- 
works and other singular diffusions, Illinois J. Math. 50 (2006) 323-383. 

[18] R. Durrett, Probability: Theory and Examples, 2nd edition, Duxbury Press, Belmont, 
California, 1996. 

[19] R. Durrett, Stochastic Calculus, CRC Press, Boca Raton, Florida, 1996. 

[20] E.B. Dynkin, Markov Processes, Vol. II, Die Grundlehren der mathematischen Wis- 
senschaften 121-122, Springer, Berlin, 1965. 

[21] S.N. Ethier and T.G. Kurtz, Markov Processes - characterization and convergence, John 
Wiley & Sons, New York, 1986. 

[22] K. Fleischmann and A. Greven, Diffusive clustering in an infinite system of hierarchically 
interacting Fisher-Wright diffusions, Probab. Theory Relat. Fields 98 (1994) 517-566. 

[23] K. Fleischmann and A. Greven, Time-space analysis of the cluster-formation in interact- 
ing diffusions, Electron. J. Probab. 1 (1996) Paper no. 6, pp. 1-46. 

[24] K. Fleischmann and J.M. Swart, Renormalization analysis of catalytic Wright-Fisher 
diffusions, Electron. J. Probab. 11 (2006) Paper no. 24, pp. 585-654. 

[25] A. Friedman, Partial Differential Equations of Parabolic Type, Prentice-Hall, Englewood 
Cliffs, 1964. 

[26] A. Greven, Renormalization and universality for multitype population models, in: In- 
teracting Stochastic Systems (eds. J.-D. Deuschel and A. Greven), Springer, Berlin, 2005, 
pp. 209-246. 

[27] F. den Hollander, Renormalization of interacting diffusions, in: Complex Stochastic Sys- 
tems (eds. O.E. Barndorff-Nielsen, D.R. Cox and C. Khippclberg), Monographs on Statis- 
tics and Applied Probability 87, Chapman & Hall, 2001, Boca Raton, pp. 219-233. 

[28] F. den Hollander and J.M. Swart, Renormalization of hierarchically interacting isotropic 
diffusions, J. Stat. Phys. 93 (1998) 243-291. 

[29] S. Kliem, Degenerate stochastic differential equations for catalytic branching networks, 
preprint. 

[30] R.G. Pinsky, Positive Harmonic Functions and Diffusion, Cambridge University Press, 
Cambridge, 1995. 

[31] S. Sawyer and J. Felsenstein, Isolation by distance in a hierarchically clustered population, 
J. Appl. Probab. 20 (1983) 1-10. 

[32] T. Shiga and A. Shimizu, Infinite-dimensional stochastic differential equations and their 
applications, J. Math. Kyoto Univ. 20 (1980) 395-416. 

[33] DW. Stroock, S.R.S. Varadhan, Multidimensional Diffusion Processes, Springer- Verlag, 
New York, 1979. 

[34] J.M. Swart, Clustering of linearly interacting diffusions and universality of their long-time 
distribution, Probab. Theory Relat. Fields 118 (2000) 574-594. 



47 



[35] J.M. Swart, Uniqueness for isotropic diffusions with a linear drift, Probab. Theory Relat. 
Fields 128 (2004) 517-524. 

[36] S.R.S. Varadhan, Probability Theory, Courant Lecture Notes 7, American Mathematical 
Society, Providence, Rhode Island, 2001. 



48 



